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We present results for the AS = 2 matrix elements which are required to study neutral kaon 
mixing in the standard model (SM) and beyond (BSM) . We also provide leading chiral order results 
for the matrix elements of the electroweak penguin operators which give the dominant AI = 3/2 
contribution to direct CP violation in K tttt decays. Our calculations were performed with 
Neuberger fermions on two sets of quenched Wilson gauge configurations at inverse lattice spacings 
of approximately 2.2 GeV and 1.5 GeV. All renormalizations were implemented non-perturbatively 
in the RI/MOM scheme, where we accounted for sub-leading operator product expansion corrections 
and discretization errors. We find ratios of non-SM to SM matrix elements which are roughly twice 
as large as in the only other dedicated lattice study of these amplitudes. On the other hand, our 
results for the electroweak penguin matrix elements are in good agreement with two recent domain- 
wall fermion calculations. As a by-product of our study, we determine the strange quark mass. Our 
main results are summarized and discussed in Sec. VII. Within our statistics, we find no evidence 
for scaling violations. 
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I. INTRODUCTION 



Flavor changing neutral currents (FCNC), and amongst them neutral kaon mixing, have been an extremely rich 
source of information about electroweak interactions and have been highly instrumental in designing the current 
standard model (SM) of particle physics. Like all FCNC in the SM, AS" = 2 processes only appear at one loop and are 
highly suppressed. Therefore, they are also a sensitive probe for new physics and provide a rich supply of constraints 
for models of physics beyond the standard model (BSM), such as supersymmetric models, left- right symmetric models, 
multi-Higgs models, etc. 

In the SM, neutral kaon mixing is dominated by VF-exchange box diagrams whose short-distance contributions 
can be described by an effective Hamiltonian in which the VF-bosons and up- type quarks have been integrated out. 
Because of the left-handed nature of the coupling of M^-bosons to quarks, only one four-quark operator arises: 

Oi = ^7^1 - 75)d''][s-V(l - 75)d'] , (1) 

where a and b are color indices. 

On the other hand, in models of new physics, where heavy right-handed particles can also experience flavor changing 
interactions, the effective AS' = 2 Hamiltonian is generically of the form: 

5 3 

n^i=' - E (/^)^' + E (/^)^' ' (2) 

where, in addition to the SM operator of Eq. (1), there are the four operators, which we give here in the supersymmetric 
basis of [1]: 

02 = r(i-75)rfi[5''(i-75)rf''] 

03 = r(l~75)rf''][s''(l-75)d"] 

04 - r(l-75)d1[5''(l+75)rf''] (3) 

05 = r(l-75)d''][s^l+75)d"] , 
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as well as the parity transformed operators [1]: 

01 = r7M(l + 75)rf"][sV(l+75)rf''] 

02 = r(l + 75)d"][s'(l + 75)rf''] (4) 

03 = r(l + 75)d''][s'(l+75)d1 • 

In Eq. (2), the Wilson coefficients Ci and Ci incorporate the short-distance physics associated with the new heavy 
degrees of freedom which may participate in extensions of the SM, as well as the SM contribution in the case of 
Ci. These coefficients can be computed perturbatively. Moreover, they depend on the flavor mixing parameters of 
the new physics model considered. These parameters can therefore be constrained by experimental measurements of 
observables sensitive to AS = 2 processes. The relevant observables include the K^-Ks mass difference, AM/^ , and 
the quantity e, which parameterizes indirect CP violation in if — > tttt decays. The short-distance contributions to 
these quantities are given by (cf. e.g. [2]): 

AMk ^ Mk, - Mks ^ Re ^^°'^y 1^°^* (5) 

T[Ks (7r7r)j=o] 2V2AMk Mk ' ^ ' 

where the subscript I denotes the isospin of the two pion state. They are measured experimentally to a precision 

better than 1% (cf. [3]). 

In addition to the short-distance contribution encoded in Eq. (5), AMk also rccciivcis "long-distance" contributions 
in the SM [4]. These arise when the quarks which participate in the box diagram involve the up quark. They can 
either be "very long-distance" , involving the exchange of pseudo-Goldstone boson pairs between two AS = 1 vertices. 
Or they can be "slightly-lcss long-distanc:e" . involving the exchange of non-Goldstonc, u. d. s hadrons. The former may 
be large, but their determination has very large uncertainties [5] . The latter have recently been studied in a model of 
large- A^c QCD and are found to be of order 10% [6] . Such contributions are not a problem in the calculation of e, since 
this quantity is dominated by the top and charm quark contributions for which the short distance approximation of 
Eq. (6) is very good. In general, assumptions have to be made about these long distance contributions, as well as 
the values of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements in the presence of new physics. Nevertheless, 
under reasonable hypotheses, interesting constraints on extensions of the SM can be found, as shown for SUSY 
extensions, for instance, in [7]. 

Regardless of the hypotheses made, to validate or constrain models of new physics with measurements of AMk 
and e, one clearly has to be able to quantify the non-perturbative, strong-interaction physics present in the matrix 
elements of the operators Oi and Oi between and states. This is where a non-perturbative technique such 
as lattice QCD enters. Now, while the operators in Eqs. (1), (3) and (4) have both parity even and odd parts, it is 
clear that only the parity even parts contribute to the matrix elements relevant for K'-^-K'^ mixing, {K'^\Oi\K^) and 
{K'^\Oi\K'^). Moreover, since the strong interaction preserves parity, we have {K'^\Oi\K'^) — {K'^\Oi\K'^) , i = 1,2,3. 
Thus, in the remainder, wc need only consider the operators of Eqs. (1) and (3). 

The transformation properties of the operators Oi under the chiral group SU{3)l x SU{3)r are interesting. It is 
commonly known that Oi belongs to the (27, 1) representation of this group, which implies that its matrix element 
between pseudoscalar states vanishes like the square of the mass of this state in the chiral limit, as can be seen from 
its vacuum saturation value given below in Eq. (7). It is straightforward to work out that O2 and O3 belong to the 
(6,6) representation, while O4 and O5 belong to the (8,8). This means that the matrix elements of these operators 
between psciidoscalar states go to a constant in the chiral limit, as indicated once again by their vacuum saturation 
values given below in Eq. (8). Thus, the matrix elements of the non-SM operators are enhanced compared to the 
chirally suppressed SM one by the large factor M^/{nis + radY' ~ 25, as obtained for current values of (mg -|- rrid) at 
2 GeV in the MS scheme. Clearly, a reliable determination of these enhanced matrix elements is needed. 

The first objective of our work is to calculate the matrix elements {K^\Oi\K^) , i = 1, • • • , 5, using lattice QCD. 
The calculations are performed using Neuberger fermions [8-11] on quenched gluon backgrounds generated with 
the standard Wilson gauge action. The gauge configurations were generated at two values of the bare coupling 
/3 = 2Nc/gQ = 6.0 and 5.85 on 18^ x 64 and 14^ x 48 lattices, respectively. With the lattice spacings determined in 
Sec. Ill, i.e. a"^ = 2.17(9) GeV at /3 = 6.0 and = 1.49(4) GeV at /3 = 5.85, these lattices have sides of size 1.6 and 
1.9 fm. This work extends to new matrix elements, larger lattices and to another value of the lattice spacing-allowing 
for an estimate of discretization errors-our earlier investigation of the SM matrix element {K^\Oi\K^) , performed on 
a 16^ X 32 lattice at (3 = 6.0 [12]. The great advantage of using Ginsparg- Wilson fermions [13] such as those proposed 
by Neuberger, is that unlike more traditional fermion formulations (e.g. staggered or Wilson fermions), they possess 
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an exact, continuum-like, chiral flavor symmetry at finite lattice spacing [14]. ^ In particular, this means that the 
operator mixing pattern is the same as in the continuum a non-negligible simplification for the operators considered 
here and that discretization errors appear only at sub-leading order, which at currently used lattice spacings is 
certainly an advantage. Moreover, to avoid potentially large perturbative corrections, we perform all renormalizations 
non-perturbatively in the RI/MOM scheme a la [16]. Preliminary results of our calculations at /3 = 6.0 were presented 
in [17, 18]. 

Since the matrix elements are quantities of dimension four, they suffer strongly from the uncertainty in the lattice 
spacing, which is of order 10% in quenched calculations [19]. It is therefore advantageous to give dimensionless 
measures of these matrix elements, by normalizing them with quantities which are similar enough that some of the 
systematic errors are expected to cancel in the ratio. For the SM contribution, determined by the operator Oi, one 
usually computes the B-parameter Bi or Bx, which measures the deviation of the matrix element from its vacuum 
saturation value. We have 

f(i^°|57,75d|0)(0|,s7^75d|i^"> fM^Fi ' 

where fj, is the renormalization scale and Mk and Fk are, respectively, the mass and decay constant of the kaon, 
which are measured precisely experimentally. Our normalization of the decay constant is such that Fk = 113 MeV 
and we take Mk = 495 MeV. 

One can also define S-parameters for the matrix elements of 02,...,5, much in the same way. Neglecting sub-leading 
chiral terms, we have: 

_ mo.ii.m _ (^°iQ,(M)|g°) 

for i = 2,- ••,5 with our phase conventions, with Ni = {|,— 2,— |} and where we have used the partial con- 
servation of the axial current, 57^75 rf, to rewrite the matrix elements of the pseudoscalar density 575 d in terms of 
the ,s and d quark masses in the second equality. However, unlike Bk which gives a dimensionless measure of the 
matrix element of Oi in terms of well measured quantities, the relation of i?2,---,5 to the corresponding matrix elements 
involves matrix elements of 575^, or equivalently light-quark masses, which are dimensionful quantities not directly 
accessible to experiment and whose determination requires a non-perturbative calculation. Thus, in providing these 
-B-parameters, we are essentially shuffling some of the systematic errors into the determination of the normalization 
factor or worse, compounding them if the correlations between the parameters and the normalization factors are not 
taken into account. This point was first emphasized in [20], where alternate suggestions for normalizing the matrix 
elements were given, all involving the mass of the K*. Since the systematic errors for vector mesons may be quite 
different from those on pseudoscalars, we propose another dimensionless measure of the matrix elements, defined 
through the ratios: 







\aP {p^\o,{^i)\p^)^ 




WkI 


ex'pt 


[f^ (pO|Oi(/x)|PO)J 


lat 



(9) 



for i = 2, • • • , 5, where M and F are the mass and "decay constant" of the lattice kaon which we denote by to 
indicated that the mass of the strange and down quarks that compose it can diff'er from their physical values. The 
factor of is present to keep the ratios from diverging in the chiral limit. Though we are not interested in this 
limit here, the strong mass dependence induced by this divergence could complicate the necessary interpolation to the 
mass of the kaon. It is worth noting that the ratios i?f ^'^(/i, M^) measure directly the ratio of BSM to SM matrix 
elements and, as such, can be used in expressions for AM^f and e beyond the SM, in which the SM contribution is 
factored out. 

Another dimensionless measure of the matrix elements which involves only quantities pertaining to pseudoscalar 
mesons are the ratios 

for i = l,---,5. From these ratios, the physical matrix elements can be obtained by considering 2Fj.Gi{iJ., M^). 
Moreover, in the chiral limit, the Gj for i = 2, • • • , 5, once multiplied by the correct Wilson coefficients, can be 



^ The same is true of domain-wall fermions [15] which are also Ginsparg- Wilson fermions in the limit of an infinite fifth dimension. 
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thought of as contributions to the coupHngs of operators in an effective, weak chiral Lagrangian. For simplicity, in 

the following wc will abusively call the Gi couplings. 

As already mentioned, O4 and O5 belong to the (8,8) representation of SU{3)l x SU(3)ii. This is the same chiral 
representation as the A/ = 3/2 parts of the electroweak penguin operators Qr and Qs, which give the dominant 
AJ = 3/2 contribution to e'. In the S'l7(3)-flavor limit, we have: 

(7r+|Q?g|i^+) = i(^0|Og,4|i^0) , (11) 

where the relative phase between the two matrix elements depends on the phase conventions chosen for pion and kaon 
states ^ and where we define: 

g?/' = i[s''74i-75)d'^][^iV(i + 75)«'-dV(i + 75)^1 

+^[s''7/.(l - 75)«1[nV(l +75)rf'] (12) 
^[s"7m(1 - 75)^1 [^^"7^(1 +75)^1 • (13) 
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In this same limit we also have 



= ^5,4 , (14) 

where i?7 g measures the deviation of {i^^lQr g \K~^) from its vacuum saturation value and here, phase choices cancel 
so that the equality of Eq. (14) is convention independent. 

Moreover, at leading order in the chiral expansion, these 7r+ matrix elements are related to the — > (7777)7=2 

matrix elements relevant for the calculation of e' in the following way: 

((7r7r)/=2|Q7,8|^^°) cx l(7r+|Q?g|i^+) + 0{p^) (15) 

where F is the pseudoscalar decay constant which, at this order, can be taken to be F^, Fk or any other value in the 
chiral regime and where the constant of proportionality depends on phase and normalization conventions. Eq. (15) 
becomes exact, of course, in the chiral limit. 

Thus, the second objective of our work is to determine the matrix elements {■k~^\Q^^^\K'^)/F, in the chiral limit. 
Ideally then, we would compute (^5^4 of Eq. (10) on the lattice for mesons P in the chiral regime and extrapolate 
the result to the chiral limit to obtain the electroweak penguin couplings G^^^'"^ of the corresponding weak, effective 
chiral Lagrangian, or equivalently, the matrix element ((7r7r)/=2|Q7,8|-f^°) in the chiral limit. However, two problems 
stand in our way. The first is that the mesons in our simulation, except for perhaps the two to three lightest ones, are 
beyond the chiral regime. The second is that the chiral behavior of 65,4 in the quenched and the Nf = 3 unquenched 
theory differ, indicating that we cannot expect quenched results for G^^^ in the chiral limit to have much to do with 
their physical values. 

Due to computational limitations, we cannot do much about the first problem at present, but we can try to address 
the second. Indeed, as discussed in Sec. VI and first argued in [18], the matrix elements {'K'^\Q^^g\K~^) can be 
normalized in such a way that the resulting quantities have the same NLO chiral logarithm in the quenched and 
Nf =3 theories. In the /S'f/(3)-flavor limit, we can consider 

D',/,'{„, M') ^ = F^ G5,4 , (16) 

where F is the decay constant common to the pions and kaons. This similarity in the chiral expansions allows for the 
possibility that the quenched and unquenched results for these quantities in the chiral limit may not be too different. 



In our conventions, the matrix elements (7r+|Qy''g |i<r+) diflfer by a sign from the ones given in [21]. 
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Because of its importance in constraining the CP violating phase of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix with the measurement of e (cf. e.g. [22, 23]), Bj^ has received much attention, both analytically and in lattice 
QCD. Analytically, it has been evaluated with two-point [24, 25] and three-point [26-28] QCD sum-rules, as well 
as in models of large A^c QCD [29-31]. The most recent calculations were performed in the chiral limit, except for 
that of [29] , where it was shown that chiral corrections represent approximatively 50% of the final number. On the 
lattice Bk also has a long history, which has been reviewed recently in [32-34]. In particular, there are at present 
two quenched benchmark calculations [35, 36] and, very recently, the first Nf = 2 + 1 unquenched calculation with 
staggered fermions, at a single value of the lattice spacing [37]. The point of our calculation is not to establish a 
new quenched benchmark for Bk, but rather to quantify K^-K^ mixing beyond the SM in a consistent framework. 
In fact, our dimensionless ratios iJ^SM q£ ^g-j^ -j^ which a number of systematic errors are expected to cancel, 
may be used in the future, in conjunction with high precision unquenched determinations of Bk, to yield improved 
calculations of K^^-K'^ mixing beyond the SM. 

Unlike the situation that we have described in the preceding paragraph for Bk, there is only one other lattice 
calculation of the full set of /S.S ~ 2 matrix elements relevant for K'^-K^ mixing [20] . It was performed with tree- level 
0(a)-improved Wilson fermions, which are not chirally symmetric and suffer for potentially large 0{asa) errors. In 
fact we find results for the BSM matrix elements which are substantially larger than those of [20], implying stronger 
constraints on BSM models than previously found. 

The situation regarding electroweak penguin matrix elements is somewhat intermediate. There are a number of 
analytical calculations [21, 38-40] as well as a handful of recent quenched lattice calculations [41-45] and a preliminary 
Nf = 2 unquenched calculation [46]. There arc also older calculations of the corresponding _B-parametcrs [47, 48] 
as well as exploratory calculations which attempt to go beyond the chiral limit [45, 46], including one [49] which 
is based on the finite- volume formalism of [50]. Yet there is no consensus on the value of these matrix elements 
even in the chiral limit. Thus, our results, obtained in a formulation of lattice QCD with full chiral symmetry, full 
0(a)-improvement, at two values of the lattice spacing and using non-perturbative renormalization, are useful for 
clarifying the current situation. 

Our paper is organized as follows. In Sec. II we briefly summarize the parameters of our simulation. In Sec. Ill 
we describe our analysis of two-point functions relevant for the determination of the matrix elements of interest. 
This includes a description of our extraction of the strange quark mass. In Sec. IV we use three-point functions to 
determine the bare i?-parameters and other quantities related to the four-quark operator matrix elements, as functions 
of our bare quark masses. Sec. V is dedicated to the non-perturbative calculation of the renormalization constants 
required to subtract the logarithmic divergences present in the bare results. In Sec. VI we discuss the quark-mass 
dependence of the results renormalized conventionally at 2GeV ^ in the RI/MOM scheme and discuss the necessary 
interpolations and extrapolations to the physical quark masses. As in the previous sections, the results obtained at 
the two values of the lattice spacing are discussed in parallel, and a direct comparison of these results is presented. 
In Sec. VII we combine the results obtained at the two values of the lattice spacing for the matrix elements, matrix 
element ratios and quark masses interpolated or extrapolated to their physical mass point. We give our final physical 
results at 2 GeV in both the RI/MOM and MS-NDR schemes, including estimates of the dominant systematic errors 
within the quenched approximation. We also compare our results to those obtained both on and off the lattice. The 
reader interested only in our final results can skip directly to this section. Finally, in Sec. VIII, we conclude with a 
few closing remarks. 

The paper further contains a number of appendices. Appendices A-D contain tables which are referred to in earlier 
sections, but which we place in appendices to improve the readability of the paper. For completeness, we provide 
in Appendix E the non-perturbative renormalization constants in the MS-NDR scheme as well as the quark-mass 
dependence of the four-quark matrix element results renormalized in the MS-NDR scheme at 2 GeV in Appendix F. 
We also give, in Appendix G, results in the MS-NDR scheme at the individual lattice spacings for the B-parametcrs 
and other matrix element ratios at the physical point. These results are used in Sec. VII to obtain our final physical 
results in that scheme. 



The renormalization scale // = 2 GeV is convenient, because it is in the perturbative domain and because 2 GeV is comparable to the 
lattice spacings that we use. The latter guarantees that we do not introduce potentially large ATj-dependent ln(a/i) logarithms in 
matching our results to continuum renormalization schemes. 
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II. NUMERICAL IMPLEMENTATION 



We describe gluons with the Wilson gauge action and quarks with the overlap action 

Q . (17) 



S = ^2 QD[mq]q = q ( 1 — ^om, j £'[0] + m, 

where niq is the bare quark mass and -D[0] is the Ncubcrgcr-Dirac operator [9] 

D[0] = ^(l + X/VxTx^ . (18) 

In Eq. (18), X = Dw — p/a is the Wilson-Dirac operator with a negative mass parameter p which separates physical 
from doubler modes. We generated 100 quenched gauge configurations each at a bare value of the gauge coupling, 
P = 6.0 and /3 = 5.85, corresponding to inverse lattice spacings of = 2.12 GeV, respectively = 1.61 GeV, as 
obtained from the Sommer scale defined by roF(ro) = 1.65 and ro = 0.5 fm [51, 52]. * On each gauge configuration, we 
invert the Ncubcrgcr-Dirac operator for a number of values of the quark mass ranging from roughly TOs/3 to 1.5tos. 
At /? = 6.0, the bare masses in lattice units are amq = 0.030,0.040,0.060,0.080,0.100 and at /3 = 5.85, they are 
aniq = 0.030,0.040,0.053,0.080,0.106,0.132. Details on our implementation of the Neuberger-Dirac operator can be 
found in [53]. All statistical errors are estimated with a single-elimination jackknife, which is propagated throughout 
our analysis. 

In our simulations, wc choose the negative mass parameter p to vary with the lattice spacing a, namely p = 1.4 
at /3 = 6.0 and p = 1.6 at (3 ~ 5.85, to optimize the locality of the Neuberger-Dirac operator [54]. At fixed lattice 
spacing, a change in p corresponds to an 0{a^) redefinition of the action and one may worry about the effect on the 
continuum limit of our procedure. However, as long as p, viewed as a function of the lattice spacing a, remains in the 
single particle sector and has a smooth limit in the range < p < 2 for a ^ 0, changing p with /3 simply represents 
a modification of the approach to the continuum limit. Moreover, for a function p which approaches its continuum 
limit value with corrections of O(a^), this modification only appears at O(a^). We require that our two values of p 
lie on such a curve. 



III. TWO-POINT FUNCTIONS, SPECTRAL QUANTITIES AND SCALE SETTING 

To obtain the vacuum saturation values of the matrix elements of interest, we compute the following two-point 
functions: 

X 

and 

41(^o) = E(^(^)^o(0))' (20) 

X 

where P = J75S is the pseudoscalar density and Aq is the time-component of the axial current A^^ = drf^^^f^s. Here 
s = (1 — aD\Q]/2p)s. The introduction of the "hatted" quark fields guarantees that the bilinear and four-quark 
operator matrix elements have the proper chiral properties and therefore renormalize as in the continuum and are free 
of 0{a) discretization errors. Moreover P and Aq are obtained from P and ^0 by interchanging the flavors d s. 
We work in the 5J7(3)-flavor limit where = = rUq and we allow the common bare quark mass am, to take on 
the values described in the preceding section. 

While the correlation function Cp'js^^{xo) is only used to normalize the three-point function involving operator Oi 

as described below, we fit the time-symmetrized correlation function Cpp{xQ) to the usual cosh form in the time 
interval 12 < xo/a < 19 at /? = 6.0 and 10 < xo/a < 14 at /3 = 5.85 for all quark mass values, obtaining the masses M 
of the corresponding pseudoscalar mesons as well as the matrix elements (0|P|P°). The initial fit time ensures that 



* As detailed in Sec. Ill, compatible values for the lattice spacing are obtained from more physical quantities, such as the leptonic decay 
constant of the kaon, Fk- 
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0.045 




M^/(4nF) 



FIG. 1: Mass-dependence at /3 = 6.0, in terms of the variable M^/(47r_F)^, of the resuhs for the decay constant in lattice units, 
aF, obtained from the two correlation functions Cpp(xo) and C^p^g^f^p_g^{xo). The solid curve is a quadratic fit to the results 

from Cpp{xo)- It is used to interpolate the results to the kaon point M'^/{AnF)^ — Mj^ /(AtvFk)'^ , shown as a vertical dotted 
line. The result of this interpolation is then used to determine the lattice spacing. 



asymptotic behavior has been reached while the final time is chosen to match the range used in three-point function 
fits, as described in Sec. IV. Because of the chiral symmetry of Neuberger fermions, the axial Ward identity can be 
used to derive the pseudoscalar decay constant F from this matrix element, without a finite renormalization. Results 
for aM and aF are summarized in Table XI for /3 = 6.0 and Table XII for (3 — 5.85. 

Because Neuberger fermions satisfy an index theorem [14, 55], one may worry about the presence of finite- volume, 

(2) 

quenched, zero mode contamination in the Cpp{xo) correlator. However, we showed in [53] that these contributions are 
not visible within our statistical accuracy at the volumes and quark masses that we consider. In [53], we compared 

(2) (2) 

directly results of fits to Cpp{xo) and C'^p_|.5)(p_5') (a;o), where S is the scalar density, and found no appreciable 
difference. The two correlators are both dominated by the pseudo-Goldstone boson state for asymptotic times xq and 
C^Ip+S)(p-S) (-^o) ^^^^ from zero-mode contributions by chirality. For completeness, this comparison is repeated here. 
Results obtained for the time-symmetrized correlation function C^p_^_g)(^p_g■^{xo), for the same time fitting ranges as 
above, are also given in Tabs. XI-XII. 

(2) 

As Tabs. XI-XII show, Cpp{xo) yields results for aM and aF which are compatible with those given by 

(2) 

^{p+s){P~s)^^^^- "^^^ compatibility is also displayed graphically in Figs. 1-2, where the results for aF from the 
two correlation functions are plotted against the corresponding values M^/(47rF)^. We conclude that the effects 
of finite-volume zero-mode contributions are not statistically significant at the volumes and quark masses that we 

(2) 

consider, at least for two-point functions at asymptotic times. For simplicity, we will consider only Cpp{xo) in the 
remainder of this paper. 

To determine the lattice spacing a, we study aF as a function of M^/(47rF)^ and interpolate aF quadratically to 
the point M^/(47rF)^ = M]^/{4ttFk)'^ , with Mk and Fk given in the Introduction. The interpolations are shown 
in Fig. 1 for /3 = 6.0 and in Fig. 2 for (3 = 5.85, where the result corresponding to aniq — 0.030 is omitted because 
MpL ~ 3.5 for this point, which makes it more susceptible to finite volume effects. With the value of the decay 
constant in lattice units obtained at the kaon point, we set fl-i = FK/{aF) and obtain a''^ = 2.17(9) GeV at /3 = 6.0 
and a~^ = 1.49(4) GeV at (3 = 5.85. These values for the lattice spacing are compatible with those derived from the 
Sommer scale in Sec. II, namely a~^ = 2.12 GeV at /3 = 6.0 and a^^ ^ 1.61 GeV at /? = 5.85, given the intrinsic 
0(10%) model-dependence associated with rg. We choose to set the scale and fix the strange quark mass with Mk 
and Fk because these quantities are closely related to the matrix elements of interest. This procedure may therefore 
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FIG. 2: Same as Fig. 1, but for /3 = 5.85. 



TABLE I: {nis + m) at 2 GeV in the RI scheme as obtained from our simulations at /3 = 6.0 and /3 = 5.85. The first error is 
statistical and the second comes from the systematic uncertainty in the determination of Zg. 

/3 (m, + m)"'(2GeV) 
6.0 120.(8)(2)MeV 
5.85 127.(6)(10)MeV 



absorb some of the quenching artefacts associated with the physics that we are studying. In any event, we will be 
considering mainly dimensionless quantities such as i3-parameters and matrix-element ratios. The value of the lattice 
spacing will then only enter logarithmically through the renormalization scale dependence of the matrix elements and 
in fixing the kaon point in the mild interpolations that we will have to perform. 

As a side product of our analysis, we determine the value of the strange quark mass. The renormalized quark 
mass is obtained from the bare quark mass in lattice units as m^^(2GeV) = (am) * a^^/Z^^(2 GeV), where the 
lattice spacing is the value determined above, and Zg^(2GeV) is the renormalization constant of the scalar density 
determined non-perturbatively in Sec. V. We assume that the meson mass depends only on the sum of the masses of 
the quarks which compose it. In Fig. 3, we show a quadratic interpolation to the kaon point of the RI/MOM sum of 
quark masses at 2 GeV, {mq^ + mg2)^^(2 GeV), as a function of M^/(47rF)^ at /3 = 6.0. Also shown in this plot are 
the results bX. (3 — 5.85, where the interpolation, performed without the point corresponding to aniq — 0.030 for the 
reasons discussed above, has been omitted for clarity. Only statistical errors are shown, and already the agreement 
is very good. The interpolation to the kaon point gives (m^ + to)^^(2 GeV), where rh is the average up and down 
quark mass. Our results for this quantity, obtained at the two values of the lattice spacing are summarized in Table 1. 
There is no evidence for scaling violations. Results in the MS scheme are given in Appendix G. The final result will 
be given in Sec. Vll. 
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0.05 0.1 0.15 0.2 0.25 

FIG. 3: Mass-dependence at /3 = 6.0, in terms of the variable M^/(47r_F)^, of the sum of the quark masses which compose the 
pseudoscalar meson of mass M. The quark masses are renormalized in the RI/MOM scheme at 2GeV. The solid curve is the 
result of a quadratic fit, plotted in the fit region. The fit is used to interpolate the results to the kaon point / {4ttF)^ = 
shown as a vertical dotted line. The dashed curve is an extension of the fit curve outside the fit range. Also 
shown are the results at /3 = 5.85, but the interpolation is omitted for clarity. Only statistical errors are shown. 



IV. THREE-POINT FUNCTIONS, BARE B-PARAMETERS AND FOUR-QUARK OPERATOR 

MATRIX ELEMENTS 

To obtain the bare matrix elements of the A^* = 2 operators, Oi, between K'^ and K'^ states, and the corresponding 
B-parameters, we construct the following ratios of three- to two-point functions: 



Es,y{P{^)Oi{0)P{y)) 



Spp(xo,yo) = g-^^ /r./ ^ J J /n^ ^^ ' ^1 (^l) 



|E.-,^(^(^)^o(0))(Ao(0)P(y)> 
Ej,,^(^(2^)O.(0)P(y)) 



Opp{xo,yo) = ^^^^^^ |-n,^^f^,n^\lr,ln^-^,t..^\ ' ' (22) 



for i = 2, • • • , 5, with the A^^'s given after Eq. (8). As indicated, these ratios become constant for asymptotically large 
initial (j/o) a-nd final (xq) times. These constants are the bare _B-parameters, i = 1, • • • ,5, whose renormalized 
analogs were defined in Eq. (8). In Eqs. (21)-(22), the subscripts on ^'j^j., * = 1, ■ • • ,5, label the sink Jf{x) and the 
source Ji{y), which we defined below Eq. (20) for the case of P{x) and Ao{x). In addition, the bare operators Oi in 
Eqs. (21)-(22) are those of Eqs. (l)-(3) with d replaced by d. 

We fit the Bpp ratios of Eqs. (21)-(22) to a constant in the symmetric time intervals given by: 12 < xo/a < 19 
and 45 < yo/a < 52 for i=l, • • • , 5, at /3 = 6.0; 10 < xo/a < 12 and 36 < yo/a < 38 for i=l and 10 < xo/a < 14 
and 34 < ya/a < 38 for i=2, • • • , 5, at P — 5.85. The points nearest the origin are are chosen so that the ratios are 
dominated by the ground state contribution. The other extremities of the ranges are picked so that time reversed 
contributions are less than 1%, even for the lightest meson (next to lightest at /3 = 5.85). The results for the bare 
_B-paranieters are given in Table XIII at /3 = 6.0 and Table XIV at /? = 5.85, as functions of bare quark mass. 

Though no evidence was found in the previous section for finite-volume, zero-mode contamination in two-point 
functions, one could worry about their effect on three-point functions. In order to verify whether these modes might 
be a problem, we considered an alternate ratio for determining Bi which, unlike Bpp{xQ,yQ) of Eq. (21), is entirely 
free of zero-mode contamination because of its chiral structure. This ratio is the one that we advocated in our first 
determination of Bj^ with Neuberger fermions [12], performed on a lattice whose four- volume is roughly three times 
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smaller than the one considered here: 

^1 , . J2s.yi^o{x)Oi{0)La{y)) a«xo«T/2<yo«T 

I Ex,? (^0 (a;)Lo (0)) (Lo (0)Lo (?/)) 

where Lq is the time-component of the left-handed weak current, = ^7^(1 — 75)5. As found in [12], this ratio 
reaches asymptotic behavior at early times. Fits of this ratio in symmetric time ranges with x^ja > 8 and yo/a< 56 
at /3 = 6.0 and with xo/a> 6 and jjn/a < 42 at /? = 5.85 gives bare values of Bi which are entirely compatible with 
those obtained from ,Bpp (xo,yo)- Wc take this as evidence that the three-point functions which appear in Eqs. (21)- 
(22) are not significantly affected by zero-mode contributions for the masses that we consider in our relatively large 
volumes. 

In computing matrix elements of the operators O2 and O3, it is also possible to find sources which entirely eliminate 
finite- volume, zero-mode contributions. In the definition of 6p'p(.xo, yo) in Eq. (22), one must replace the pscudoscalar 
source and sink in the numerator hy R = S + P, whore S is the scalar density, S = ds. This eliminates zero-mode 
contributions in the three-point function by chirality. In the denominator, the two point function {P{x)P{Q)} must 
be replaced by {R{x)L{0)) , and (P(O)P(j/)) by (Z(0)i?(j/)), where L = S — P. Both these two-point functions are 
free from zero modes. The resulting ratios, B'^^{xo, yo) are thus also free from zero-mode contributions. As one can 
easily convince oneself, such a manipulation is not possible in correlation functions relevant for computing the matrix 
elements of the operators O4 and O5. There, zero- mode contributions cannot be eliminated by simply using chirality. 
Nevertheless, in [18], we have checked that sources with different chiral properties, and thus different zero-mode 
contributions, give identical results. In any case, we do not pursue this issue further here and assume that the absence 
of zero-mode contributions observed in pseudoscalar two-point functions and in the three-point function of Eq. (21) 
is evidence that these contributions are not a problem. And we take as our final results for the P-parameters those 
given by the ratios of Eqs. (21)-(22), because the use of identical sources for the different operator matrix elements 
yields cancellations of some statistical fluctuations in the ratio of these matrix elements. 

(2) 

From these B-parameters and our determinations of M and F from the two-point function Cpp, we then derive a 
number of other quantities of interest. The building blocks for these quantities are the bare matrix elements of our 
flve operators Oj, i = 1, • • • , 5, which we obtain as follows: 



1 fi 

(pO|Oi|P°) = -^M^F^B, (24) 
(P°|0,|P") = -NmP\P°)\^Bi . (25) 



In Eq. (24), Za is the normalization constant of the local axial current, A^, defined after Eq. (20). It is there because 
the vacuum saturation value of the matrix element of Oi, used in the definition of Bi, is obtained in terms of the 
matrix clement of Aq, while in constructing the matrix element in Eq. (24), we have used the partially conserved 
axial current, through the axial Ward identity. Z-^ could be omitted here and reintroduced when we renormalize our 
results. However, for our simulation parameters it leads to a large correction (see Table II). We have thus chosen to 
introduce it. because it allows for direct comparison of quantities containing (P''|Oi|P'^) of Eq. (24) to those that are 
obtained using a vacuum saturation calculated with matrix elements of the local axial current ^o- 

In fact, we have performed such a comparison using (O|^o|-P*') obtained from both Cp\ and C^^p, together with 

(2) 

Cpp. The bare matrix elements determined using the three different estimates of its vacuum saturation value are 

entirely compatible. However, the matrix element's statistical error does depend quite strongly on the procedure 

(2) 

used to obtain it from Bi. The matrix element obtained using a vacuum saturation value from Cp^^^ is the one with 

(2) 

the smallest statistical error. It is followed by the one obtained which makes use of C\^p and finally by the one of 

Eq. (24), which makes use of only Cpp. Nonetheless, the best measure of the physical matrix element {K^\Oi\K'^) 
is given by -^M|-P|-Pi, which does not necessitate reconstructing (P°|Oi|P°) for our lattice mesons, P". Moreover, 

obtaining (P°|Oi|pO) from Bi and a fit to c'pl^ gives results for the BSM ratios of Eq. (9) which have significantly 
larger statistical errors than those obtained using Eq. (24). The same is true of the matrix element obtained from 

Bi and a fit to C)^^p. In fact, both these determinations of BSM ratios have very comparable errors. We thus take 

(P^IOilP") as obtained through Eq. (24) and avoid discussing fits of the two-point functions Cp\^ and C^^p which 
have already been detailed in [53]. 

In addition, as discussed in the Introduction, we have chosen not to give the values of the matrix elements themselves, 

which have mass dimension four and are therefore very sensitive to the intrinsic uncertainty associated with the 
determination of the lattice spacing in quenched calculations. Rather we choose to present two dimensionless measures 
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TABLE II: Results for the renormalization constant Za of the local axial current at /? = 6.0 and /3 = 5.85. 

/3 Za 
6.0 1.5548(11) 
5.85 1.4434(18) 



of these matrix elements, namely the BSM ratios of Eq. (9) and the couplings of Eq. (10), the latter enabhng a 
determination of the matrix element itself for a given meson mass, using the corresponding value of the leptonic decay 
constant. 

Thus, from the matrix elements obtained using Eqs. (24)-(25) and the values of M and F from fits to C);^ (Tabs. XI- 
XII), we straightforwardly construct the bare equivalents of the BSM ratios defined in Eq. (9). We summarize our 
results for these ratios as a function of bare quark mass at /? = 6.0 in Table XV and at /3 = 5.85 in Table XVI. 

We also construct the bare equivalents of the couplings Gi of Eq. (10), from the matrix elements as given by 
Eqs. (24)-(25), and the pseudoscalar decay constant obtained from Cpp and the axial Ward identity and given, 
again, in Tabs. XI-XII. Our results for these contributions to the bare couplings as a function of bare quark mass are 
summarized in Table XVII for /3 = 6.0 and Table XVIII for /3 = 5.85. 

Finally, in view of obtaining the chiral limit values of the matrix element of the electroweak penguin operators 
relevant for direct CP violation vnK ^ -k-k decays, we determine the quantities D^^^^ defined in Eq. (16). Our results 
for these quantities in lattice units are summarized in Table XIX for /3 = 6.0 and Table XX for /? = 5.85. 



V. NON-PERTURBATIVE RENORMALIZATION 



The bare AS' = 2 matrix elements of operators Oj, i = 1,...,5, have logarithmic divergences which must be 
subtracted. Because Neuberger fermions exhibit, at finite lattice spacing, a chiral flavor symmetry analogous to that 
of contimmm QCD. the renormalization goes through as in the contimmm. In particular, Oi, the standard model 
operator renormalizes multiplicatively, while O2 and O3 mix under renormalization, as do O4 and O5. These mixing 
properties are also compatible with the chiral properties of the operators described in the Introduction. 

To avoid potentially large perturbative errors, we compute all renormalization constants non-perturbatively in the 
RI/MOM scheme a la [16]. Thus, we fix gluon configurations to Landau gauge and numerically compute the relevant 
amputated forward quark Green functions with legs of four- momentum p, Aq^ {'mq,p, go), where nig is the mass of the 
quark. The renormalization constants are then determined by requiring that the corresponding renormalized vertex 
functions have their tree-level values. We therefore define the ratios: 

7?RICrr, n\ - v2 Tt {Ay {nig , p, go)Vvf 

{mg,p,go) = Z^—j- — , (26) 

TT[{Ao^{mq,p,go)Voi\ 

where i — j = 1 for Oi, i,j G {2, 3} to determine the mixing of O2 and O3 and i, j G {4, 5} for the mixing pair O4, O5. 
In Eq. (26), the Vo are normalized projectors onto the spin-color structure of the tree- level operator O [56]. Because 
of the chiral symmetry of Neuberger fermions, the parity even and odd parts of the operators renormalizc in the same 
way. We therefore consider the parity odd parts, as CPS symmetry then guarantees that off-diagonal mixing terms 
vanish exactly. 

Chiral symmetry also guarantees that the product Za x Tr {Ay (mg,p)'Py} is the quark wavefunction renormal- 
ization squared, Zg. Since we can obtain Za very cleanly from the axial Ward identity [53], this method allows us 
to circumvent a more direct calculation of Zg, a procedure which induces large discretization errors. We choose Ay 
instead of A^ because the former has a slightly cleaner signal. In Table II, we summarize the results for Za obtained 
from the axial Ward identity. 

The momentum combinations that we consider in constructing the ratios of Eq. (26) are the following. We allow 
Pi, 2, 3 X {L/2tt) to take on values from to 2, and po x (T/27r), values from to 4. For larger momenta, we restrict 
Pi, 2, 3 X {L/2-k) to take on values from 2 to 3 and pq x [T/^-k), values from 7 to 10. Moreover, we use the lattice-fermion 
definition of momentum, p^ = {1/a?) X];i=o (^^^ '^Pfj.Y- These choices reduce 0(4)-breaking discretization errors so that 
the TZ^^{mg,p, go) are to a good approximation functions of p^. 

At large (p^ > 1.5 GeV^), an operator product expansion (OPE) in l/p"^ can be performed on the Green functions 
which appear in the definition of TZf^ {nig ,p,go)- In this expansion, the leading term, proportional to the unit operator, 
gives the desired renormalization constants in the RI/MOM scheme. Perturbative mass effects appear at sub-leading 
orders, suppressed by powers of p"^, as do possible non-perturbative contributions. On the lattice, of course, 0(4) 
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symmetry is broken and there will be discretization errors. However, as already noted, for our choice of momenta, we 
expect the Ti-fj^ {mqjp, go) to be functions of to good approximation. Therefore, the leading terms in the OPEs of 
these ratios at finite niq, including discretization error terms, are of the form (see also [12]): 

n^{mq,p,9o) = ■■■ + ^S'(m,^^9.,p') ^ ^RGi(^^^ + B^\g„p-) [apf + • • • , (27) 

where the ellipses stand for higher-order terms in the OPE and sub- leading discretization errors. Here U^{p'^) 
describes the running of the renormalization constants in the RI/MOM scheme, Zjj^{p'^) = Ujj^{p'^)Z^j'^ , where the 
superscript "RGI" (renormalization group invariant) indicates that the corresponding renormalization constant is the 
constant which translates lattice results into scale- and scheme-independent quantities. ^ The p^ dependence in A^^^ 
and B^j"^ is specified to accommodate the expected logarithmic corrections to the power behavior in p^. Note that 

is independent of quark mass, because this dependence only appears at higher orders in the OPE. 
To determine whether power corrections and discretization errors are relevant in extracting the renormalization 
constants from the TZ^^{mq,p, go), we consider the "RGI" ratios: 

^S^'(m„p,5o) = [U''%\p'')Kj{mq,p,go) (28) 

obtained from the TZ^^ {mq,p, go) by dividing out the running. In the absence of discretization errors and sub- leading 
perturbative and OPE terms, these "RGI" ratios should be constant functions of p^. In Fig. 4, we plot TZfj^^{mq,p, go) 
for i = j = 1 as functions of for our smallest and largest quarks masses, arUq = 0.03 and 0.10 at /3 = 6.0. We do 
the same in Fig. 5 for TZ^'^^{niq,p, go) with i, j G 2,3, and in Fig. 6 for the mixing of operators O4 and O5. Here 
the running is implemented at two loops [57, 58], with Nf = and Aqcd = 0.238(19) GeV from [59]. A number of 
observations are in order: 

• to a very good approximation, TZ^^^ {mq,p, go) are functions of p^; 

• the TZfj^^{mq,p, go) are not generically constant functions of p^ at values of where one might consider ex- 
tracting the renormalization constants; when they are, this behavior seems to be the result of a turnover instead 
of genuine p^-independence; 

• the ^Z^j^^{mq,p, go) do in some cases exhibit a significant mass dependence, but that dependence tends to vanish 
when p^ increases as predicted by the OPE; 

• discretization errors are clearly visible but remain moderate and there is no evidence for contributions of order 
higher than (ap)^. 

Our conclusion is that power corrections and discretization errors may give significant contributions and they should 
be accounted for in extracting the renormalization constants, as already briefly discussed in [12]. This is even more 
true at /3 = 5.85 where discretization errors enter at lower physical momenta. So, instead of performing, in some cases 
delicate, chiral extrapolations of TZ^^{mq,p, go) to mq = ^ and taking Z^^{go,p^) = TZ^^{'mq,p, go) for p^ ~ 4 GeV^ 
as is usually done [16, 56], we fit the TZ^^ {mq,p, go) at finite ruq to OPE forms such as the one given in Eq. (27). 
This procedure not only eliminates the mass effects through their p^ dependence, but also other power-suppressed 
contributions as well as discretization errors. In performing these fits, we neglect possible logarithmic corrections in 
p^ in the coefficients ^^j'^' and -B-J-*. They are at least subdominant to the leading corrections and sub-subdominant 
compared to the perturbative behavior that we want to isolate. ^ Moreover, they are not known analytically and 
they would be difficult to identify in the data. We also neglect discretization errors in (am)^ since they should be 
significantly smaller than either the (ap)^ discretization errors or the 1/p^ corrections in the range of p^ considered. 

We obtain our central values for the renormalization constants by fitting to all points p^ such that 3 GeV^ < P^ < 
14GeV^ at /3 = 6.0 and 3GeV^ < P^ < 8.5 GeV^ at /3 = 5.85. The reason for this rather large range of p^ is that we 



^ To specify our normalization conventions for RGI quantities, we give as an example the factor [f/'^']!! at one loop. It is (2/3oas/47r)4/2/5o 

where, of course, /80 = H — 2Arj/3 is the coefHcient of the one loop beta functions. 
® We have performed the calculation at five quark masses which are not that close to the chiral limit, and the full mass dependence of 

these functions, especially in the quenched approximation, is not known. 

l/p** corrections are also formally sub-subdominant. However our results often clearly indicate the presence of such terms. 
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FIG. 4: Plot of the renormalization ratio TZi'i{mq,p, go) for the standard model operator Oi as a function of four-momentum 
squared, p^, for our most massive (am, = 0.10, circle) and our lightest (am, = 0.03, cross) quarks. The solid curve for each 
data set illustrates the combined fit to the parameterization of Eq. (27), with the terms indicated in Table III, and in range 
fitted, as described in the text. The dotted curve is an extension of this fit beyond the fit range. The solid horizontal line is 
the resulting value of Zf{^^{go). 
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FIG. 5: Same as Fig. 4, but for the operator pair {O2, O3}. Here all eight ratios are fit simultaneously. 



want to make sure that all OPE and discretization terms, which make a significant contribution to TZfj^{mq,p, go) in 
regions where we can hope to extract reliable values of the renormalization constants, are accurately measured with 
the data. Moreover, down to 3GeV^, the two-loop running remains under control. For TZ^f{mq,p, go), we consider 
even lower momenta, down to = 1.5 GeV^. The reason for this is the presence of a large power correction whose 
coefficient can only be determined at these low values of p^. And because the running of Oi is mild, the two- loop 
evolution of this operator remains sensible even down to these low values of the momentum. The terms kept in the 
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FIG. 6: Same as Fig. 4, but for the operator pair {O4, O5}. Here all eight ratios are fit simultaneously. 



TABLE III: Terms retained in the fits of the renormalization ratios TZf^{mq,p, go) to the OPE form of Eq. (27), in the range 
3 GeV^ < < 14 GeV^ at /? = 6.0 and 3 GeV^ < < 8.5 GeV^ at = 5.85. The discretization term in (ap)^ is treated as mass 
independent, as discussed in the text. At /3 = 6.0, [am'^'" , am';^'''') = (0.030,0.100) and (am';^'"' ,am';;^°'^) = (0.040,0.132) at 
f3 — 5.85. TZ^^ is the ratio from which the renormalization constant of the scalar density is obtained. It is fit to an OPE form 
starting from 2GeV. 
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fits for the various renormalization constants are summarized in Table III. 

We initially individually fit TZ^^{mq,p, go) obtained for our lightest or next to lightest and heaviest quark masses 
((aTO™",oTO^°^) = (0.030,0.100) at /3 = 6.0 and (0.040,0.132) at /3 = 5.85). We find that, despite a factor of 
over 3 between these two masses, the resulting Zjj'^^{go), as well as the discretization error coefficients Bij{go), are 
independent of quark mass within statistical errors, giving us confidence in the consistency of our OPE description. 
Therefore, for our final results we perform a combined fit to the TZfj^{mq,p, go) for the two masses, constraining 

^ij'^^i9o) and Bij{go) to be mass- independent, but allowing, of course, A[^^ to be different for the two masses. The 
resulting fits at /? = 6.0 are shown in Figs. 4-6, together with the horizontal line describing the behavior of the 
perturbative term alone, as reconstructed from the fit. The distance between this line and the data points quantifies 
the impact of OPE and discretization corrections to the leading perturbative term. Similar plots are obtained at 
f3 — 5.85, though the stability of the fits is reduced by the presence of discretization errors at even lower values of p^. 
The resulting renormalization constants, Z^^{2 GeV), in the RI/MOM scheme at 2 GeV are summarized in Table XXI 
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for /3 = 6.0 and Table XXII for (3 ~ 5.85. Since the MS-NDR scheme is also often used, for completeness we give the 
renormalization constants in that scheme in the Appendix E in Tabs. XXXVIII-XXXIX. 

To renormalize the _B-parameters, we further need to renormalize the matrix elements' vacuum saturation values. 
For 02,- -.5, these are given by matrix elements of the pseudoscalar density, so that the renormalization constant needed 
is Zp. This constant is actually difficult to obtain directly in the RI/MOM scheme of [16], because the relevant vertex 
receives chirally enhanced power corrections [16]. Fortunately, the chiral symmetry of Neuberger fermions guarantees 
that Zp = Zs, where Zs is the renormalization constant for the scalar density, which is free of Goldstone pole 
contaminations. Thus, we consider a ratio TZ^^{mq,p, go), which is analogous to the ratios TZf^^{mq,p, go) defined in 
Eq. (26), with the operators Oi and Oj replaced by the scalar density and the squares replaced by first powers in 
the numerator. As above, we study the momentum dependence of this ratio and fit it to an OPE expression of the 
form given in Eq. (27), in the interval 2 GeV^ < < 14GeV^ at /? = 6.0 and 2 GeV^ < < 8.5 GeV^ at /? = 5.85. 
We find that a good description of the data is obtained by retaining the discretization term proportional to (ap)^ 
together with 1/p^ and 1/p^ OPE terms (see Table III). The data and the fit at /? = 6.0 are shown in Fig. 7, together 
with the horizontal line describing the behavior of the perturbative term alone, as reconstructed from the fit. One 
difference here is that the running of the operator is implemented at four loops [60-62], with Nf = and the same 
value of Aqcd as above. The result in the RI/MOM scheme at 2 GeV is given in Table XXIII. The value at the same 
scale in the MS scheme is given in Appendix E in Table XL. These constants are also required for determining the 
renormalized quark masses of Tabs. I-LV, from the corresponding bare quark masses. 

Though it is straightforward to obtain the renormalization constants for the i3-parameters from those for the 
operators in Tabs. XXI-XXII and for Zp in Table XXIII, we nevertheless give our results for these constants in the 
RI/MOM scheme at /3 = 6.0 in Table XXIV and at /3 = 5.85 in Table XXV, because they may be of use and because 
we can account for the correlations in the ratios of Zij to Zp. Results in the MS-NDR scheme are summarized in 
Tabs. XLI-XLII in Appendix E. 

To estimate the systematic errors associated with our procedure for determining renormalization constants, we 
consider fits to a higher momentum range, where power corrections in 1/p^ will play a less important role. We choose 
the range 7.5 GeV^ < < 14GeV^ at /3 = 6.0 and 5GeV^ < < 8.5 GeV^ at /? = 5.85, where momenta belong to 
the high momentum cone described below Eq. (26). We also eliminate the power corrections in 1/p^ to which our 
results are not sensitive in this higher range of momenta. The list of terms kept is summarized in Table IV and the 
resulting renormalization constants are given in Tabs. XXVI-XXVII in the RI/MOM scheme. The resulting values 
for Zs are given in Table XXIII and the results for the renormalization constants of the i?-parameters are summarized 
in Tabs. XXVIII-XXIX. 

As the renormalization constants obtained from the two fitting ranges indicate, our description of the RI/MOM 
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TABLE IV: Same as Table III, but for the restricted range 7.5 GeV^ < < 14GeV^ at /3 = 6.0 and 5GeV^ < < 8.5 GeV^ 
aX 13 = 5.85. 
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Green functions is very stable at /3 = 6.0, giving us confidence that we have isolated correctly the various OPE and 

discretization terms which contribute to these functions. However it is clear that for a high precision calculation 
of matrix elements, a non-perturbative renormalization procedure in which corrections in powers of 1/p^ are absent 
and in whic;li pcrtm-bation theory is only used at much higher scales would be highly desirable, especially on coarser 
lattices. Proposals in that direction have been made very recently in [63, 64] . 



VI. MASS DEPENDENCE AND PHYSICAL MATRIX ELEMENTS 



Having determined the renormalization constants, wc can now renormalizc the bare matrix elements and B- 
parameters for our various meson masses. Renormalizing these quantities before studying their behavior in mass 
is useful because it allows for comparisons of this behavior with results obtained in other simulations, performed with 
either different actions or different parameters. In our case, wc undertake such a comparison between the results 
obtained on our finer and coarser lattices, at /? = 6.0 and /? = 5.85, respectively. Moreover, this renormalization is 
also important if one is to perform a combined chiral and continuum extrapolation using finite lattice spacing chiral 
perturbation theory. 

In Tabs. XXX-XXXl we present the results obtained at /? = 6.0 and f3 — 5.85, respectively, for the B-parameters in 
the Rl/MOM scheme at 2 GeV, as functions of bare quark mass. The corresponding results for the renormalized BSM 
ratios, Rf^^, are given in Table XXXII for /3 = 6.0 and Table XXXIII for f3 = 5.85. In Table XXXIV for /3 = 6.0 
and Table XXXV for (3 = 5.85, we summarize our results for the couplings Gi of Eq. (10), again as a function of 
light-quark mass. Finally, Tabs. XXXVI-XXXVII display the results that we obtain for D^^^ in the RI/MOM scheme 
at /? = 6.0 and f3 = 5.85, respectively. For completeness, we present all of these results renormalized in the MS-NDR 
scheme at 2 GeV in Appendix F, in Tabs. XLVII-LIV. 

In order to finalize our calculation of the S-parameters, BSM ratios Rf^^ and couplings Gi relevant for K^-K^ 
mixing beyond the standard model, we have to interpolate our renormalized results to the kaon mass. To perform 
the interpolation wc fit our data, with the exclusion of the single point for amq = 0.030 at /? = 5.85 which we leave 
out for the reasons explained in Sec. HI, to simple polynomial forms in M"^ / {At: F)"^ . As we did for setting the scale 
in Sec. Ill, we fix the kaon point to be M^/{4:TtF)^ = M]^/ [AttFkY . More complicated chiral perturbation theory 
fits could be tried, but their applicability in our mass range is questionable and for the simple interpolations that we 
have to perform, any reasonable description of the data is sufficient. 

We begin with the interpolation of the S-parameters. We find that and are consistent with linear behavior 
in M^/ [A-kFY while B2, B3 and B4 are well fit by parabolas in our range of masses. We show these fits for illustrative 
purposes in Fig. 8, but do not reproduce the fit parameters here since they are not of physical interest. Our results 
for the S-parameters at the kaon mass in the RI/MOM scheme at 2 GeV are summarized in Table V at /3 = 6.0 and 
5.85. 

We now turn to the BSM ratios Rf^^. Here we find that R^^^ and i?^^ are weU fit by a straight Hue while R^^^ 
and Rf^^^ require quadratic terms in A/p /{AttF)^ . The results of the interpolations to the kaon point are reproduced 
in Table V and are illustrated in Fig. 9. 
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FIG. 8: Mass-dependence, in terms of the variable M^/(47r_F)^, of the B-parameters Bi, i = 1, • ■ ■ , 5, in the RI/MOM scheme 
at 2 GeV. The solid curves are the results of the fits described in the text, and are plotted in the fit region. The fits are used 
to interpolate the results to the kaon point M^/{4ttF)^ — M]^ /(AttFk)'^ , shown as a vertical dotted line. The dashed curves 
are an extension of the fit curves outside the fit range. 
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FIG. 9: Same as Fig. 8, but for the BSM ratios Rf 



2,- ■■,5. 



We further interpolate the couphngs Gi. As Fig. 10 confirms, Gs is well fit by a constant, G4 by a straight line, 
while Gi, G2 and G3 require quadratic terms in M^/(47ri^)^. Again, we do not reproduce the fit parameters here. 
We only note that in the chiral limit Gi is consistent with zero, which is the expected behavior. The results for the 
couplings Gi at the kaon mass in the RI/MOM scheme at 2 GeV are given in Table V. 

From these results for Bi, Gi and /if^'^ at the kaon mass in the RI/MOM scheme, it is straightforward to obtain 
the corresponding results in the MS-NDR scheme at 2 GeV. This is done in Appendix G and the results are presented 
there in Table LVI. 

Finally, we turn to the most delicate part of this discussion: the chiral extrapolation of the electroweak penguin 
matrix elements. As stated in the Introduction, the quantities D^^^ of Eq. (16) have the same chiral logarithm at 
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FIG. 10: Same as Fig. 8, but for the "couplings" d, i = 1, • • ■ ,5. Also shown, for indicative purposes, is an extrapolation of 
Gi to the chiral limit, which is denoted by a dashed vertical line at A/^/(47r_F)^ — 0. 



TABLE V: Results at the kaon mass in the RI/MOM scheme at 2GeV for fJ — 6.0 and /3 = 5.85. The first error is statistical 
and the second originates from the systematic uncertainty in the renormalization constants. 





/3 = 6.0 






f3 = 5.85 




i B, 






B, 


G. 




1 0.563(47)(7) 


28.8(25)(3) 


1 


0.534(36) (4) 


27.3(18)(0) 


1 


2 0.865(72)(9) 


-464.(96)(11) 


-16.1(30)(6) 


0.898(89)(44) 


-430. (69) (40) 


-15.8(25)(14) 


3 1.41(10)(6) 


151.(30)(4) 


5.24(93) (9) 


1.53(10) (40) 


146.(17)(13) 


5.37(64)(51) 


4 0.938(48) (13) 


601.(87)(0) 


20.7(27) (2) 


0.904(73) (106) 


514.(69)(24) 


18.8(27)(7) 


5 0.616(51)(14) 


128.(11)(1) 


4.57(60)(1) 


0.56(14)(1) 


108.(27)(17) 


3.9(11)(6) 



NLO in the quenched and Nf = 3 theories, while the corresponding matrix elements or parameters G5.4 do not. 
Indeed, at NLO in the Nf ^ 3 theory [65], 

and [66] 



-+I0?;8l^^> I 1-477^ In 3AM,.. +W ) ■ (30) 



In the quenched theory, we have at NLO [67], 
and [68] 

(.-|0?5|A-) ^ A<?5- (1 - ^ In J^^^^ + 0(p-,) , (32, 
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FIG. 11: Mass-dependence, in terms of the variable M^/{4-KFf, of the ratios D^^g^ in the RI/MOM scheme at 2GeV. The 
solid curves are the results of the fits described in the text, and are plotted in the fit region. The fits are used to interpolate 

shown as a vertical dotted line. The dashed curves are an extension 



the results to the kaon point Af^/(47rf )^ = A/^/(47rFif 

30 shown, 

is denoted by a dashed vertical line at M^/{-iTTF)^ 



of the fit curves outside the fit range. Also shown, for indicative purposes, are extrapolations of D^^g to the chiral limit, which 

0. 



where it is understood that the low energy constants pp, ■ ■ ■ may differ in the quenched and Nf — 3 theories. Thus, 

3 /2 

for Z?7 g , we find in both cases, the following chiral behavior: 

<^<-(l-^...^^j;^.0(.^,). ,33, 

where again, the low energy constants are distinct in the two theories. However, the fact that the functional behavior 
is the same and the fact that the coefficient of the logarithm is small suggests that D'^^^'^ may be comparable in the 
two theories, especially since our lattice has been tuned to reproduce results at the kaon mass and the extrapolation 

3/2 

should be mild. This mildness of the mass behavior is actually seen in our results for Dy^ renormalized in the 
Rl/MOM scheme at 2 GeV, which are plotted as functions of AP/iinF)^ in Fig. 11. 

3/2 

Instead of attempting a hazardous chiral extrapolation of Dj^ , we choose to interpolate these quantities to the 

kaon point. Because Fk is used to fix the lattice cutoff, we take D^^^p,, M|-) = F]^ x G54(/i, M|.), with G5a{^, ^k) 
obtained above in the RI/MOM scheme at 2 GeV. We then use the fact that any reasonable chiral extrapolation (linear, 
quadratic, log form of Eq. (33), log form plus a quadratic term in M^/(47rF)^, . . .) gives values which are compatible 

3/2 . . 3/2 

with the results at the kaon point. This is shown in Fig. 11 where Dg is extrapolated linearly and Z?g quadratically 
for illustration. Thus, the systematic uncertainty due to the chiral extrapolation should be smaller than our statistical 
errors and we have approximately D^^^'^ ~ D'^^^{M'^). Now, for determining from D'^^^'^ the chiral limit value of 

the K ^ TTTT matrix element required to calculate e', we can use Eq. (15), taking ((7r7r)/=2|Q7,8|-f^")x ^ Px^^[s ^^ 
with F^ = 87(5) MeV [65]. Alternatively, since we are working to leading chiral order, we can choose to remain at the 

kaon point and consider ((7r7r)7=2|Q7,8|-?^'')x ^ ^kD^^^^M]^). For our final value of the matrix elements, we take an 
average of the two central values and use their spread as an estimate of the chiral extrapolation error. This leads to a 
±15% systematic error on ((7r7r)/=2|Q7.8|^")x- Table VI summarizes all of these results in the RI/MOM scheme at 
2 GeV for /3 = 6.0 and 5.85. The results in the MS-NDR scheme are reproduced in Table LVII of Appendix G. 
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TABLE VI: Results, at (3 = 6.0 and 5.85, for the ratios -D^^f of Eq. (16) interpolated to the kaon point and our best estimate 
for ((7r7r)/=2|Q7,8|A'") in the chiral limit written as a matrix element according to Eq. (15), to avoid ambiguities due 

to normalization conventions for the two pion state. All results are given in the RI/MOM scheme at 2 GeV. Please see text for 
a discussion of how the central value and systematic error on the matrix elements is obtained. 





D^'^iMl) [GeV^] 


Dl'\Ml) [GeV^] 




[GeV^] 


(^+\qI'''\k+) 


[GeV^] 


F 


F 


6.0 
5.85 


1.63(14)(1) 
1.38(34)(21) 


7.7(11)(0) 
6.6(9)(3) 


0.163(14)(1)(21) 
0.138(34)(21)(18) 


0.77(11)(0)(10) 
0.66(9) (3) (9) 



TABLE VIL Final results at the kaon mass in the RI/MOM scheme at 2 GeV. The first error is statistical and the second is 
systematic and is described in the text. Also shown for comparison are the results of [20], obtained using quenched, tree-level, 
0(o)-improved Wilson fermions. 



i B, 






Ri 




This work [20] 


This work 


[20] 


This work 


[20] 


1 0.563(47) (30) 0.69(21) 


28.8(25) (16) 


37.(9) 


1 


1 


2 0.865(72)(35) 0.70(9) 


-464. (96) (35) - 


-0.24(3) X 10^ 


-16.1(30)(6) 


-6.4(18) 


3 1.41(10)(13) 1.1(1) 


151. (30) (7) 


74.(9) 


5.24(93)(16) 


2.0(6) 


4 0.938(48)(36) 1.1(1) 


601.(87)(87) 


4.4(4) X 10^ 


20.7(27)(18) 


12.(3) 


5 0.616(51)(59) 0.77(11) 


128.(11)(20) 


1.0(2) X 10^ 


4.57(60) (68) 


2.8(8) 



VII. FINAL PHYSICAL RESULTS AND DISCUSSION 

As seen the preceding Section, results at p = 6.0 and 5.85 are in excellent agreement. We could therefore assume 
that we are in the scaling regime and fit the dependence of our results to a constant. However, as discussed briefly 
in Sec. V, the determination of the renormalization constants is difficult at /? = 5.85. We therefore choose to take the 
results obtained on our finer lattice as our final central values and use the results at /? = 5.85 to estimate possible 
discretization errors. The lattice spacing squared at /3 = 6.0 is essentially half way between that at /3 = 5.85 and the 
continuum limit. Assuming that a linear extrapolation of the results in a? to the continuum limit gives a worst case 
scenario, wc estimate our symmetric discretization errors to be the absolute value of the difference of the central values 
of the results obtained at the two lattice spacings. We then add this error in quadrature with the error coming from 
the systematic uncertainty in the renormalization constants. The procedure assumes that 0{a?) errors are dominant 
over higher order discretization errors and that their coefficient is the same at the two values of /3, as discussed in 
Sec. II. 

Let us first consider the BSM matrix elements. We give our final results for the B-parameters, couplings Gi and 
BSM ratios Ri in Table VII, in the RI/MOM scheme at 2 GeV. In this table we also provide the results derived from 
the only other calculation of these matrix elements [20]. The latter were obtained in the quenched approximation, 
using tree-level, 0(a)-improvcd Wilson fermions. The first point that we wish to make is that our result for Bk = Bi 
is in excellent agreement with the recent world averages of lattice results given in [32-34]. The second point is that our 
BSM matrix elements are enhanced at Mk by a factor of five to twenty compared to the chirally suppressed standard 
model one, as shown by the values that wc obtain for the BSM ratios Ri. The third point is that this enhancement 
is significantly larger than the one found in [20]. This difference can only partially be explained by our smaller value 
of B\, which accounts for about a 23% enhancement in our results for Ri over those of [20]. The remainder of the 
enhancement must therefore come from elsewhere. It does not come from the _B-parameters Bi for i = 2, • • • , 5, 
which agree within errors. Since the _B-parameters agree much better than do the matrix elements themselves, we 
must already disagree on the vacuum saturation values, that is on the matrix element of the pseudoscalar density, 
(0|s75C?lif°). Since the axial ward identity requires that (0|s75d|i^°) = \/2FKM]^/{ms + m), we must in fact find a 
smaller (m^ -|- m). Using the results of Tabs. I-LV and the method described above for obtaining final results on the 
matrix elements, we find, for the quark masses: 

{rus + to)^^(2 GeV) = 120.(8)(7) MeV ^ (m^ + mf^{2 GcV) = 102.(7)(6) MeV , (34) 
which is in excellent agreement with the earlier Neuberger fermion calculations of [53, 69]. Our result for this sum 
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TABLE VIII: Final results at the kaon mass in the NDR scheme at 2GeV. 
systematic and is described in the text. 



The first error is statistical and the second is 



Ri 



0.571(48)(30) 29.2(26)(16) 1 
0.679(56)(27) -512.(106)(39) -17.5(32)(7) 
1.055(77)(98) 159.(31)(7) 5.44(97)(15) 
0.810(41)(31) 730(105)(106) x 10^ 24.7(32)(22) 
0.562(39)(46) 165.(15)(25) 5.76(73)(78) 



TABLE IX: Final results for {{7r7r)i=2\Q7,8\K'') in the chiral limit written as a K'^ — > tt"*" matrix element according to Eq. (15), 
so as to avoid ambiguities due to normalization conventions for the two pion state. Results are given in the RI/MOM and 
MS-NDR schemes at 2 GeV. Please see text for a discussion of how the central value and systematic error on the matrix 
elements is obtained. 



scheme 


(7r+|Q^/^|if+> 


[GeV^] 


(7r+|Q^/ = |K+> 


[GeV^] 


F 


F 


RI/MOM 
NDR 


0.163(14)(33) 
0.211(20)(42) 


0.77(11)(15) 
0.93(13) (18) 



of quark masses is almost 28% lower than the one found by the authors of [20], using tree-level, 0(a)-improved 
Wilson fcrmions at the couplings which they use to obtain the BSM matrix elements. Indeed, in [70], they obtain 

{mg+fh)^'^^{2 GeV) = 130(2) (18) MeV with the same gauge configurations as in [20]. Since this mass appears squared 
in the relation of the BSM matrix elements to their respective S-parameters, this difference accounts for the large 
enhancement of our BSM matrix elements over those of [20] and, combined with our lower value of Bi , for the factor 
of roughly two enhancement of our BSM ratios Ri. Note that the result of Eq. (34) is compatible with the rigorous 
lower bounds calculated in [71, 72]. 

For completeness, we determine the value of the strange quark mass. Using nis/rh = 24.4 ± 1.5 as obtained from 
chiral perturbation theory [73] , the results in Eq. (34) correspond to 

mfi(2 GeV) = 116. (8) (7) MeV ^ mf^{2 GeV) = 98. (6) (6) MeV , (35) 

where the uncertainty in the chiral perturbation theory mass ratio does not contribute significantly to the error. 
Our result is fully compatible with the quenched, continuum limit, benchmark result of the ALPHA collaboration, 
m^^(2GeV) = 97.(4) MeV, obtained using Fk, as we do, to set the lattice spacing [74]. We therefore believe that the 
results of [20] for the BSM matrix elements suffer from large systematic errors which are absent in our calculation. 

Because the MS-NDR scheme is frequently used, we also provide our results for the B-parameters, couplings Gi 
and BSM ratios in this scheme at 2 GeV in Table VIII. 

We turn now to the matrix elements of the electroweak penguin operators, ((7r7r)/=2|(57,8|^°)- We compute them in 
the chiral limit, as discussed in Sec. VI. We obtain final results, as for the BSM matrix elements and quark masses, by 
combining our results at /3 = 6.0 and 5.85 in the manner discussed above. These results are summarized in Table IX, 
both in the RI/MOM and MS-NDR schemes at 2 GeV. 

These matrix elements have already been obtained, both in quenched lattice QCD and with a variety of analytical 
techniques. In Table X, we compare our results with those obtained by others in the NDR scheme at 2 GeV, also in the 
chiral limit. The results quoted for CP-PACS [42] were obtained by multiplying their results for ((7777)7=2 IQ7, 8 1-?^°) 
at 1.3 GeV in the MS-NDR scheme by — \/2/3(i^;^ -l-i^/f )/(2F^), to match our definition of the matrix elements. They 
were then run up from 1.3 to 2 GeV using Nf = two loop running with the parameters discussed in Appendix E. 

Those of RBC [43] were obtained from the expression {'k+\Q^(^\K+)/F^ = {F^ + FK){h\f^^^ J p) x a"^ with 

/ = 137(10) MeV [75]. These bare results were then renormalized in the RI/MOM scheme with the rcnormalization 
constants given in the paper at 2.13 GeV, matched onto the MS-NDR scheme and finally run down to 2 GeV with our 
two-loop running formulae. Finally, the results from SPQcdR [45] are obtained by extrapolating, to the chiral limit, 
quenched results for — > tttt matrix elements computed for a variety of initial and final state kinematics, using NLO 
chiral perturbation theory. 
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TABLE X: Comparison of our results for the matrix elements of the electroweak penguin operators in the chiral limit, with 
those obtained in other lattice calculations (above the double line) and with a variety of analytical techniques (below the double 
line). The "??" indicate that we have not estimated systematic errors in converting the results to our conventions. Moreover, 
the factors of a/2/3 in front of the SPQcdR results [45] are required, we believe, to translate their results to our normalization. " 
Results are given in the MS-NDR scheme at 2 GeV. 



Ref. 


Action 




[GeV^] 




[GeV^] 


F 


F 


This work 

CP-PACS'Ol [42] 

RBC'02 [43] 
SPQcdR'04 [45] 


Neuberger 
Domain- Wall 
Domain- Wall 
Wilson 


0.211(20)(42) 

0.220(6) (??) 
0.255(12)(??) 
.^273 X 0.16(3) 


0.93(13)(18) 
0.92(3)(??) 
1.02(4)(??) 
^y2/3 X 0.82(15) 



Bijnens et al '01 [39] 

Cirigliano ot al '02 [40] 
Friot et al '04 [76] 
Knecht et al '01 [21] 
Narison '00 [38] 



0.24(3) 
0.22(5) 
0.12(2) 
0.11(3) 
0.21(5) 



1.2(7) 
1.50(27) 
2.00(36) 
3.5(1.1) 
1.40(35) 



"'We thank C.-J. David Lin for correspondence on this subject. 



Agreement of our results with lattice calculations performed using domain- wall fermions is excellent for both matrix 

elements. Only the Wilson fcrmion results of [45] differ significantly. This may be due to the absence of chiral symmetry 
in the Wilson discretization which significantly complicates the rcnormalization procedure and allows potentially large 
0{a) discretization errors. In regards to the results obtained using analytical techniques, agreement is found with the 
results of [38, 39] and to a lesser extent [40]. The authors of [21, 76], however, find results for the matrix clement of 

3 /2 . .3/2 

which are a factor of two lower than ours, and for the matrix element of Qg , a factor of two to three higher. 
It should be noted that a quenched investigation of electroweak penguin matrix elements with modified overlap 
fermions, and perturbative renormalization, has been performed in [44]. However, the functional form for the chiral 
extrapolation of the matrix elements is incorrect and at the kaon mass, the results obtained at two different lattice 
spacdngs do not appear to be c;ompatible. In addition, preliminary results from a quenched calculation using HYP 
staggered fermions were presented in [77]. Finally, SPQcdR has also attempted to determine the ((7r7r)/=2]Q7,8|^'^°) 
matrix elements at NLO in chiral perturbation theory by fitting the required low energy constants to quenched lattice 
results for a variety of initial and final state kinematics [45] . 



VIII. CONCLUSION 



We have presented results for the full set of AS = 2 matrix elements which are required to study K^-K^ mixing 

in the standard model and possible extensions. We have also presented results, at leading chiral order, for the 

3 /2 

matrix elements of the electroweak penguin operators ^ , which give the dominant Al = 3/2 contribution to e'. 
As a by-product of our calculation, the strange quark mass was obtained. Quarks were simulated with Neuberger 
fermions, which possess an exact chiral flavor symmetry at finite lattice spacing. This not only greatly simplifies the 
complicated mixing pattern under renormalization of the operators that we consider, but also guarantees that our 
results are free of the leading 0{a) discretization errors. To avoid potentially large perturbative errors, we implemented 
all renormalizations non-perturbatively in the RI/MOM scheme, modifying the techniques of [16] to accommodate 
non-perturbative power corrections and high-momentum discretization errors. For most renormalization constants, 
we find that these effects are important. The calculations were performed on two sets of quenched configurations, 
generated with the Wilson gauge action on an 18^ x 64 lattice at /3 = 6.0 and on a 14^ x 48 lattice at /3 = 5.85. The 
two lattice spacings allowed us to investigate scaling violations. Within our statistics, we find no evidence for such 
violations and use the results on our coarser lattice to estimate possible discretization errors. 

Our main conclusion is that the non-SM, AS = 2 matrix elements are significantly larger than found in the only 
other dedicated lattice study of these amplitudes [20] . In tracing the source of this difference, we found that we already 
disagree on the much simpler matrix element of the pseiidoscalar density between a kaon state and the vacuum, which 
is the building block for the vacuum saturation values of the BSM A^* = 2 amplitudes. Through the axial Ward 
identity, the matrix element of the pseudoscalar density is related to the sum of the strange and down quark masses, 
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which wc find to be roughly 30% smaller than the value obtained in [70], with the same tree- level improved Wilson 
fermion action and gauge configurations as used in [20] . Since our result for this sum of masses is in agreement with 
the continuum limit, benchmark result of [74], we are convinced that the stronger enhancement of non-SM AS = 2 
matrix elements that we observe is correct. 

3 /2 

Regarding Bk and the matrix elements of the electroweak penguin operators Q^'g , we find good agreement with 
recent quenched lattice calculations, in particular those which make use of domain- wall fcrmions [35, 42, 43, 78], 
which have an approximate chiral symmetry at finite lattice spacing. Agreement with analytical results depends on 
the particular calculation. 

Beyond this, our investigation validates the use of Ginsparg- Wilson fermions for the calculation of weak matrix 
elements. The simplified operator mixing and non-perturbative 0(a) improvement guaranteed by the full chiral flavor 
symmetry of these fermions has proven essential for obtaining reliable results in this delicate calculation. We hope 
to have convinced the reader that these benefits outweigh the numerical overhead that accompanies the use of these 
fermions. Of course, our results suffer from the shortcomings of the quenched approximation, though we have tried to 
minimize their impact by considering dimensionless measures of the various matrix elements. Nevertheless, because 
of the numerous advantages already discussed, we believe that the way to remove this approximation is by studying 
Ginsparg- Wilson quarks on unquenched backgrounds generated with numerically cheaper fermion formulations, such 
as 0(a)-improved Wilson fermions. Work in that direction is beginning. 
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APPENDIX A: RESULTS FOR THE PSEUDOSCALAR MESON MASSES AND DECAY CONSTANTS 

AT /3 = 6.0 AND 5.85 

In this appendix, wc have regrouped the results for the pseudoscalar meson masses and decay constants, in lattice 
units, which are discussed in Sec. III. Please refer to that section for explanations. 



TABLE XI: Pseudoscalar meson mass and decay constant, in lattice units, as a function of bare quark mass at /3 = 6.0. 





{PP) 


({P + S){P 


-s)} 


aniq aM 


aF (M/AnFf 


aM aF 


(M/A-KFf 



0.030 0.215(5) 0.0514(16) 0.111(9) 

0.040 0.244(4) 0.0529(14) 0.135(9) 

0.060 0.295(3) 0.0559(13) 0.176(9) 

0.080 0.339(3) 0.0590(12) 0.209(9) 

0.100 0.379(2) 0.0620(12) 0.236(9) 



0.212(8) 0.0504(14) 0.112(8) 

0.242(6) 0.0521(14) 0.137(7) 

0.295(4) 0.0556(13) 0.178(7) 

0.340(4) 0.0592(13) 0.209(7) 

0.381(3) 0.0628(13) 0.234(8) 



TABLE XII: Pseudoscalar meson mass and decay constant as a function of bare quark mass at /3 = 5.85 in lattice units. 





{PP) 


{{P + SKP- 


-s)) 


arUq aM 


aF {M/4nFf 


aM aF 


{M/4nFf 



0.030 0.254(7) 0.0727(22) 0.077(6) 

0.040 0.287(6) 0.0738(19) 0.096(6) 

0.053 0.326(5) 0.0755(18) 0.118(6) 

0.080 0.395(4) 0.0795(16) 0.156(6) 

0.106 0.454(3) 0.0836(15) 0.187(6) 

0.132 0.507(3) 0.0876(14) 0.212(6) 



0.235(14) 0.0713(26) 0.069(12) 

0.275(11) 0.0721(21) 0.092(10) 

0.320(9) 0.0740(19) 0.118(8) 

0.395(6) 0.0789(18) 0.159(6) 

0.456(5) 0.0838(17) 0.188(6) 

0.512(4) 0.0887(16) 0.211(6) 
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APPENDIX B: RESULTS FOR THE BARE B-PARAMETERS, BSM RATIOS, "COUPLINGS" d AND 
RATIOS Dj{^ AS A FUNCTION OF BARE QUARK MASS AT /3 = 6.0 AND 5.85 

In this appendix, we have regrouped the results for the bare B-parameters, BSM ratios, "coupHngs" Gi and ratios 

3/2 

£>7 g as a function of bare quark mass at /3 = 6.0 and 5.85. These results are discussed in Sec. IV. Please refer to 
that section for explanations. 



TABLE XIII: Bare B-parameters as a function of bare quark mass at /? = 6.0. 

arriq Bi B2 -B3 B4 

0.030 0.632(63) 0.697(54) 0.803(55) 0.802(50) 0.793(50) 

0.040 0.651(51) 0.690(44) 0.797(46) 0.855(44) 0.861(45) 

0.060 0.682(40) 0.693(34) 0.801(35) 0.923(35) 0.962(39) 

0.080 0.705(35) 0.703(28) 0.810(29) 0.966(28) 1.040(33) 

0.100 0.727(32) 0.716(24) 0.822(26) 0.996(24) 1.105(30) 



TABLE XIV: Bare B-parameters as a function of bare quark mass at /3 = 5.85. 



aruq B\ B2 Bs B4 B5 

0.030 0.588(73) 0.823(49) 0.999(64) 0.640(42) 0.635(44) 

0.040 0.617(55) 0.794(43) 0.948(53) 0.721(39) 0.725(42) 

0.053 0.634(43) 0.773(38) 0.910(45) 0.794(36) 0.811(40) 

0.080 0.661(34) 0.755(30) 0.872(35) 0.887(33) 0.937(37) 

0.106 0.686(30) 0.753(26) 0.861(29) 0.940(31) 1.025(37) 

0.132 0.709(28) 0.759(23) 0.862(25) 0.976(30) 1.094(37) 



TABLE XV: Bare BSM ratios Rf^^ as a function of bare quark mass at /3 = 6.0. 

am, i?fSM ^BSM ^BSM ^BS=^ 



0.030 -19.6(32) 4.51(72) 27.0(39) 8.9(13) 

0.040 -16.6(21) 3.84(47) 24.7(29) 8.29(95) 

0.060 -13.4(12) 3.10(26) 21.4(17) 7.45(61) 

0.080 -11.59(79) 2.67(17) 19.1(12) 6.86(43) 

0.100 -10.37(60) 2.38(13) 17.32(88) 6.40(34) 



TABLE XVI: Bare BSM ratios Rf^^ as a function of bare quark mass at ^ = 5.85. 



„^ dBSM oBSM nBSM oBSM 

uTTlq 1x2 -n-3 1x4 -£X5 

0.030 -20.7(33) 5.01(78) 19.3(33) 6.4(11) 

0.040 -16.9(21) 4.05(49) 18.5(24) 6.19(82) 

0.053 -14.5(14) 3.42(33) 17.9(17) 6.09(60) 

0.080 -11.67(79) 2.70(19) 16.5(11) 5.80(39) 

0.106 -10.05(55) 2.30(13) 15.06(82) 5.47(30) 

0.132 -8.98(41) 2.039(97) 13.85(66) 5.18(25) 
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TABLE XVII: Bare ratios Gi as a function of bare quark mass at /3 = 6.0. 

arUq Gi G2 G3 G4 Gs 

0.030 12.2(19) -262.(47) 60.(10) 362.(55) 119.(18) 

0.040 15.3(19) -229.(32) 52.9(70) 340.(42) 114.(14) 

0.060 20.9(19) -194.(19) 44.8(42) 309.(27) 107.5(95) 

0.080 25.7(20) -173.(14) 39.9(30) 285.(19) 102.5(70) 

0.100 29.9(21) -160.(11) 36.6(24) 266.(15) 98.5(56) 



TABLE XVIII: Bare ratios Gi as a function of bare quark mass at /3 = 5.85. 

arUq Gi G2 G3 G4 G5 

0.030 9.2(15) -298.(47) 72.(11) 278.(41) 92.(14) 

0.040 11.9(15) -257.(30) 61.4(71) 280.(31) 94.(11) 

0.053 15.1(14) -226.(20) 53.3(46) 279.(23) 94.9(84) 

0.080 20.9(16) -190.(12) 43.8(28) 267.(16) 94.2(62) 

0.106 25.9(17) -169.4(95) 38.7(21) 254.(13) 92.2(53) 

0.132 30.4(18) -156.4(80) 35.5(17) 241.(11) 90.2(46) 



TABLE XIX: Bare ratios D^^g in lattice units as a function of bare quark mass at /3 = 6.0 . 

2 T^''>/2 2 n3/2 

aniq a a DJ 

0.030 0.96(11) 0.316(34) 

0.040 0.951(87) 0.319(28) 

0.060 0.966(61) 0.336(21) 

0.080 0.992(45) 0.356(16) 

0.100 1.023(36) 0.378(14) 



TABLE XX: Bare ratios D^^^ in lattice units as a function of bare quark mass at /3 = 5.85. 



2 7-v3/2 2 r->3/2 

arUq a Dg' a Dj' 

0.030 1.47(17) 0.487(55) 

0.040 1.53(13) 0.511(43) 

0.053 1.59(10) 0.541(35) 

0.080 1.690(74) 0.595(27) 

0.106 1.773(65) 0.644(26) 

0.132 1.853(62) 0.693(26) 
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APPENDIX C: TABLES OF RENORMALIZATION CONSTANTS IN THE RI/MOM SCHEME AT /? = 6.0 

AND 5.85 

In this appendix, wc have regrouped the tables of renormahzation constants discussed in Sec. V. Please refer to 
that section for explanations. 



TABLE XXL Mixing matrix for the operators Oi in the RI/MOM scheme at 2 GeV for /3 = 6.0. 

i/j 1 2 3 4 5 

1 2.125(20) 

2 1.851(55) 0.022(9) 

3 -0.141(7) 2.760(33) 

4 1.771(43) -0.468(21) 

5 -0.168(6) 2.541(18) 



TABLE XXIL Mixing matrix for the operators in the RI/MOM scheme at 2 GeV for /3 = 5.85. 



1 1.750(24) 

2 1.872(205) -0.133(21) 

3 -0.217(12) 2.220(53) 

4 1.939(184) -0.358(100) 

5 -0.224(10) 2.174(35) 



TABLE XXIII: Results for the renormalization constants of the local scalar and pseudoscalar densities in the RI/MOM scheme 
at 2 GeV, Z^\2GeV) = ^f^(2GeV), for /3 = 6.0 and /3 = 5.85. The results are obtained from fits to an extended range of 
as well as to a more restricted range, as described in the text. 

f3 Extended p'^-range Restricted p^-range 
6.0 1.229(14) 1.221(13) 
5.85 1.285(34) 1.194(30) 



TABLE XXIV: Mixing matrix for the B-parameters Bi in the RI/MOM scheme at 2 GeV and at /3 = 6.0. 

i/j 1 2 3 4 5 

1 0.879(8) 

2 1.225(26) -0.074(30) 

3 0.019(1) 1.828(45) 

4 1.173(16) -0.930(54) 

5 -0.037(2) 1.683(36) 
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TABLE XXV: Mixing matrix for the S-parameters Bi in the RI/MOM scheme at 2 GeV and at /3 = 5.85. 



i/j 


1 


2 3 


4 


5 


1 


0.840(12) 











2 





1.133(98) 0.403(53) 








3 





0.026(2) 1.34(10) 








4 








1.174(82) 


-0.650(201) 


5 








-0.045(3) 


1.316(87) 



TABLE XXVI: Mixing matrix for the operators Oi in the RI/MOM scheme at 2 GeV and at /3 = 6.0. It is obtained from a 
fit to a limited range and is used to determine the systematic error associated with our non-perturbative renormalization 
procedure. 



i/j 


1 


2 3 


4 


5 


1 


2.151(18) 











2 





1.806(40) 0.011(19) 








3 





-0.146(12) 2.788(33) 








4 








1.774(42) 


-0.468(20) 


5 








-0.176(10) 


2.550(17) 



TABLE XX\'ll: Sauu- as Tabk' XX\4. l)ut for J = o.cSr). 



i/j 



1 1.742(27) 

2 1.689(116) -0.185(73) 

3 -0.230(37) 2.251(48) 

4 1.869(129) -0.416(64) 

5 -0.234(31) 2.179(33) 



TABLE XXVIII: Same as Table XXVI, but for the B-parameters Bi at P = 6.0. 



i/j 


1 


2 


3 


4 


5 


1 


0.890(8) 














2 





1.212(18) 


-0.036(63) 








3 





0.020(2) 


1.870(50) 








4 











1.190(17) 


-0.941(56) 


5 











-0.039(2) 


1.711(38) 



TABLE XXIX: Same as Table XXVI, but for the B-parameters Si at /3 = 5.85. 

i/j 1 2 3 4 5 

1 0.836(13) 

2 1.184(55) 0.647(245) 

3 0.032(5) 1.579(109) 

4 1.311(63) -0.875(174) 

5 -0.055(8) 1.529(96) 



29 



APPENDIX D: RESULTS FOR THE B-PARAMETERS, BSM RATIOS, "COUPLINGS" d AND RATIOS 
L>7 8 , RENORMALIZED IN THE RI/MOM SCHEME, AND AS A FUNCTION OF BARE QUARK MASS 

AT /? = 6.0 AND 5.85 

3 /2 

In this appendix, we have regrouped the results for the B-parameters, BSM ratios, "coupUngs" Gi and ratios D^'^ , 
renormaUzed in the RI/MOM scheme at 2 GeV, and as a function of bare quark mass at /? = 6.0 and 5.85. These 
results are discussed in Sec. VI. Please refer to that section for explanations. 



TABLE XXX: B-parameters in the RI/MOM scheme at 2 GeV as a function of bare quark mass at /3 = 6.0. 

aruq Bi B2 -B3 B4 B5 

0.030 0.556(56) 0.869(74) 1.42(11) 0.911(59) 0.588(51) 

0.040 0.572(46) 0.861(62) 1.405(90) 0.971(53) 0.654(51) 

0.060 0.599(37) 0.865(48) 1.413(73) 1.046(43) 0.761(50) 

0.080 0.620(32) 0.876(42) 1.429(63) 1.094(35) 0.852(50) 

0.100 0.639(29) 0.892(38) 1.449(58) 1.127(30) 0.932(51) 



TABLE XXXI: B-parameters in the RI/MOM scheme at 2 GeV as a function of bare quark mass at /3 = 5.85. 

arriq B\ B2 B3 B4 Bs 

0.030 0.494(62) 0.96(11) 1.67(12) 0.723(67) 0.42(11) 

0.040 0.518(48) 0.925(96) 1.59(11) 0.814(71) 0.49(13) 

0.053 0.533(38) 0.900(89) 1.534(99) 0.895(75) 0.55(14) 

0.080 0.555(31) 0.879(82) 1.476(92) 0.999(82) 0.66(15) 

0.106 0.576(28) 0.876(80) 1.461(90) 1.058(86) 0.74(16) 

0.132 0.595(26) 0.883(80) 1.464(89) 1.096(88) 0.81(17) 



TABLE XXXII: BSM ratios Rf^'^ in the RI/MOM scheme at 2GcV as a function of bare quark mass at /3 = 6.0. 



am, rT"-" RF""' fif"'" 

0.030 -17.3(28) 5.65(89) 21.8(31) 4.70(70) 

0.040 -14.7(18) 4.81(58) 19.9(22) 4.47(53) 

0.060 -11.9(10) 3.89(32) 17.3(13) 4.19(37) 

0.080 -10.27(71) 3.35(22) 15.39(89) 3.99(30) 

0.100 -9.19(54) 2.99(16) 13.93(68) 3.84(26) 



TABLE XXXIII: BSM ratios Ri in the RI/MOM scheme at 2 GeV as a function of bare quark mass at /3 = 5.85. 



am, R^^ R^^ J^P^ Rf^^ 

0.030 -22.7(43) 7.9(13) 20.5(38) 4.0(13) 

0.040 -18.6(31) 6.42(80) 19.7(31) 3.9(12) 

0.053 -16.0(23) 5.44(54) 19.0(26) 3.9(11) 

0.080 -12.8(16) 4.31(32) 17.5(22) 3.84(98) 

0.106 -11.0(13) 3.68(23) 16.0(19) 3.72(89) 

0.132 -9.9(12) 3.27(18) 14.7(17) 3.60(81) 
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TABLE XXXIV: Gi ratios in the RI/MOM scheme at 2 GeV as a function of bare quark mass at /3 = 6.0. 

amq Gi G2 G3 G4 G5 

0.030 26.0(41) -493.(88) 161.(27) 621.(92) 134.(20) 

0.040 32.6(41) -431.(61) 141.(19) 584.(70) 131.(16) 

0.060 44.4(43) -365.(37) 119.(11) 530.(44) 128.(12) 

0.080 54.6(45) -326.(27) 106.3(82) 488.(31) 127.(10) 

0.100 63.6(47) -300.(22) 97.6(65) 455.(23) 125.5(89) 



TABLE XXXV: Gi ratios in the RI/MOM scheme at 2 GeV as a function of bare quark mass at /? = 5.85. 



Gi G2 G3 G4 



0.030 16.0(27) -574.(112) 200.(31) 519.(90) 101.(31) 

0.040 20.9(26) -494.(80) 170.(20) 522.(76) 104.(30) 

0.053 26.4(26) -435.(62) 148.(13) 519.(67) 107.(29) 

0.080 36.6(29) -365.(46) 122.5(84) 497.(60) 109.(27) 

0.106 45.3(32) -326.(40) 108.6(67) 471.(55) 110.(26) 

0.132 53.2(34) -301.(36) 99.7(56) 448.(52) 110.(24) 



TABLE XXXVI: D^^g ratios in the RI/MOM scheme at 2 GeV and in units of GeV^ as a function of bare quark mass at /3 = 6.0. 



n3/2 n3/2 
^8 ^7 



0.030 7.8(13) 1.67(27) 

0.040 7.7(11) 1.73(25) 

0.060 7.80(93) 1.89(23) 

0.080 8.01(82) 2.08(23) 

0.100 8.25(78) 2.27(23) 



TABLE XXXVII: D^^g ratios in the RI/MOM scheme at 2 GeV and in units of GeV^ as a function of bare quark mass at 
/? = 5.85. 



„3/2 ^3/2 



0.030 6.10(101) 1.18(35) 

0.040 6.31(93) 1.25(36) 

0.053 6.57(90) 1.35(37) 

0.080 6.98(92) 1.53(39) 

0.106 7.31(96) 1.70(41) 

0.132 7.63(100) 1.87(43) 
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APPENDIX E: NON-PERTURBATIVE RESULTS FOR THE RENORMALIZATION CONSTANTS IN 

THE MS-NDR SCHEME 

Because they may be of use in future calculations, we give non-perturbativc results for the renormalization constants 
in the MS-NDR scheme. These are obtained from the RGI renormalization constants obtained with the fits described 
in Sec. V, through the relation Zi^^^{2GeY) = Ull^^{2GeV)Zfp\ Here, U^l^^ip"^) describes the running of the 
renormalization constants in the MS-NDR scheme which, again, we implement at two loops [58] with Nj = and AcjcD 
from [59]. The matching of the renormalization constants for the pseudoscalar density Zp = Zg is implemented at 
four loops [60-62], with Nf = and the same value of Aqcd- Our results for all of these renormalization constants are 
given in Tabs. XXXVIII-XLVI, including our estimate of the systematic error associated with the analysis performed 
in Sec. V. 



TABLE XXXVIII: Mixing matrix tor the opcrattirs (), in tlic AlS-NDl? scliomc at 2 GeV for /3 = 6.0. 

i/j 1 2 3 4 5 

1 2.155(21) 

2 2.060(61) 0.007(9) 

3 -0.081(8) 2.834(33) 

4 2.152(53) -0.346(22) 

5 -0.209(8) 2.494(18) 



TABLE XXXIX: Mixing matrix for the operators Oi in the MS-NDR scheme at 2 GeV for /? = 5.85. 

i/j 1 2 3 4 5 

1 1.774(24) 

2 2.080(228) -0.152(22) 

3 -0.179(13) 2.280(54) 

4 2.356(224) -0.227(101) 

5 -0.277(12) 2.128(34) 



TABLE XL: Results for the renormalization constants of the local scalar and pseudoscalar densities in the MS scheme at 2 GeV, 
Z^^(2GeV) = ^1*^(2 GeV), for /3 = 6.0 and f3 = 5.85. The results are obtained from fits to an extended range of as well as 
to a more restricted range, as described in the text. 

f3 Extended p^-range Restricted p^-range 
6.0 1.457(16) 1.447(16) 
5.85 1.524(40) 1.416(36) 



TABLE XLI: Mixing matrix for the B-parameters Bi in the MS-NDR scheme at 2 GeV and at /? = 6.0. 



1 0.891(9) 

2 0.970(21) -0.018(22) 

3 0.008(7) 1.335(33) 

4 1.014(14) -0.489(37) 

-o.o;i;i(i) i.i7i(25) 
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TABLE XLII: Mixing matrix for the B-parameters Bi in the MS-NDR scheme at 2 GeV and at /? = 5.85. 

i/j 1 2 3 4 5 

1 0.852(12) 

2 0.896(78) 0.328(38) 

3 0.015(1) 0.982(73) 

4 1.015(71) -0.293(141) 

5 -0.040(3) 0.916(61) 



TABLE XLIII: Mixing matrix for the operators Oi in the MS-NDR scheme at 2 GeV and at /3 = 6.0. It is obtained from a 
fit to a Umited range and are used to determine the systematic error associated with our non-perturbative renormaUzation 
procedure. 



i/j 


1 


2 


3 


4 


5 


1 


2.181(19) 














2 





2.010(44) 


-0.004(19) 








3 





-0.085(14) 


2.862(34) 








4 











2.156(51) 


-0.345(21) 


5 











-0.219(12) 


2.502(17) 



TABLE XLIV: Same as Table XLIII but for /3 = 5.85. 

i/j 1 2 3 4 5 

1 1.767(27) 

2 1.874(128) -0.204(76) 

3 -0.194(41) 2.311(49) 

4 2.271(157) -0.288(70) 

5 -0.288(37) 2.132(33) 



TABLE XLV; Same as Table XLIII but for the B-parameters Bi at P = 6.0. 

i/j 1 2 3 4 5 

1 0.902(8) 

2 0.959(14) 0.010(46) 

3 0.008(1) 1.366(36) 

4 1.029(15) -0.494(39) 

5 -0.035(2) 1.194(26) 



TABLE XLVI: Same as Table XLIII but for the B-parameters Bi at = 5.85. 

i/j 1 2 3 4 5 

1 0.848(13) 

2 0.935(43) 0.508(180) 

3 0.019(4) 1.153(79) 

4 1.133(55) -0.432(125) 

5 -0.048(7) 1.064(67) 
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APPENDIX F: RESULTS FOR THE B-PARAMETERS, BSM RATIOS, "COUPLINGS" d AND RATIOS 
Z)7 8 , RENORMALIZED IN THE MS-NDR SCHEME, AND AS A FUNCTION OF BARE QUARK MASS 

AT /? = 6.0 AND 5.85 

Because MS-NDR is a popular scheme, we give in this appendix results for the B-parameters, BSM ratios, "cou- 
plings" Gi and ratios D^^, in that scheme at 2 GeV, and as a function of bare quark mass at /3 = 6.0 and 5.85. While 
wc do not use these results in our paper, they should be useful for comparisons with future lattice calculations done 
in the MS-NDR scheme. The results are obtained by combining the bare results of Sec. IV, which are tabulated in 
Appendix B, and the renormalization constants obtained in Appendix E, using the extended range. 



TABLE XLVII: S-parameters in the MS-NDR scheme at 2 GeV as a function of bare quark mass at /9 = 6.0. 

arriq B\ B2 B3 B4 B^ 

0.030 0.564(57) 0.682(58) 1.060(80) 0.787(51) 0.539(41) 

0.040 0.580(46) 0.676(49) 1.051(67) 0.839(46) 0.593(40) 

0.060 0.608(37) 0.679(38) 1.057(54) 0.904(37) 0.679(38) 

0.080 0.629(33) 0.688(33) 1.069(47) 0.945(30) 0.749(37) 

0.100 0.648(30) 0.701(30) 1.084(43) 0.974(26) 0.810(37) 



TABLE XLVIII: B-parameters in the MS-NDR scheme at 2 GeV as a function of bare quark mass at /3 = 5.85. 

aruq Bi B2 B3 B4 B5 

0.030 0.501(63) 0.752(83) 1.251(90) 0.624(58) 0.395(83) 

0.040 0.526(48) 0.726(76) 1.192(79) 0.703(61) 0.453(91) 

0.053 0.540(39) 0.707(70) 1.147(73) 0.773(65) 0.511(99) 

0.080 0.563(31) 0.690(65) 1.104(67) 0.863(71) 0.60(11) 

0.106 0.584(28) 0.688(63) 1.093(66) 0.914(74) 0.66(12) 

0.132 0.604(26) 0.693(63) 1.095(65) 0.946(76) 0.72(12) 



TABLE XLIX: BSM ratios Rf^^ in the MS-NDR scheme at 2 GeV as a function of bare quark mass at /3 = 6.0. 



amq flf"^ 

0.030 -18.9(31) 5.87(92) 26.1(37) 5.97(86) 

0.040 -16.0(20) 4.99(60) 23.9(26) 5.63(63) 

0.060 -12.9(11) 4.03(34) 20.7(16) 5.18(42) 

0.080 -11.18(77) 3.47(23) 18.4(11) 4.87(33) 

0.100 -10.00(59) 3.10(17) 16.68(81) 4.63(28) 
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TABLE L: BSM ratios Rf^^ in the MS-NDR scheme at 2 GeV as a function of bare quark mass at /3 = 5.85. 



jBSM 



0.030 -24.7(47) 8.2(13) 24.6(46) 5.2(14) 

0.040 -20.3(33) 6.66(83) 23.6(37) 5.1(13) 

0.053 -17.4(25) 5.64(56) 22.8(32) 5.0(11) 

0.080 -14.0(18) 4.47(33) 21.0(26) 4.9(10) 

0.106 -12.0(15) 3.82(24) 19.1(23) 4.64(91) 

0.132 -10.7(13) 3.39(18) 17.6(21) 4.44(83) 



TABLE LL Gi ratios in the MS-NDR scheme at 2 GeV as a function of bare quark mass at /3 = 6.0. 



G2 G3 G4 G5 



0.030 26.3(42) -544.(97) 169.(28) 754.(111) 172.(25) 

0.040 33.0(42) -476.(67) 148.(20) 709.(84) 167.(20) 

0.060 45.1(44) -403.(41) 125.(12) 643.(53) 161.(14) 

0.080 55.3(46) -360.(30) 111.8(86) 593.(37) 157.(11) 

0.100 64.5(47) -332.(24) 102.6(68) 553.(28) 153.3(98) 



TABLE LIL G, ratios in the MS-NDR scheme at 2 GeV as a function of bare quark mass at (3 = 5.85. 



arUq Gi G2 Gs G4 Gs 

0.030 16.2(28) -633.(124) 210.(33) 630.(109) 133.(33) 

0.040 21.2(27) -545.(89) 179.(21) 634.(92) 136.(31) 

0.053 26.8(26) -480.(68) 156.(14) 631.(82) 139.(30) 

0.080 37.1(29) -402.(52) 128.7(89) 604.(73) 140.(28) 

0.106 45.9(32) -359.(44) 114.1(70) 572.(67) 139.(26) 

0.132 54.0(34) -332.(40) 104.8(59) 544.(63) 137.(25) 



TABLE LIIL D^^g ratios in the MS-NDR scheme at 2 GeV and in units of GeV^ as a function of bare quark mass at /3 = 6.0. 



amq 


r,3/2 

^8 




0.030 


9.4(15) 


2.15(34) 


0.040 


9.3(14) 


2.20(31) 


0.060 


9.5(11) 


2.37(29) 


0.080 


9.7(10) 


2.57(27) 



0.100 10.02(94) 2.78(28) 



TABLE LIV: D^ g ratios in the MS-NDR scheme at 2 GeV and in units of GeV^ as a function of bare quark mass at /3 = 5.85. 

am, J Dl ~ Dj '' 

0.030 7.4(12) 1.56(38) 

0.040 7.7(11) 1.65(38) 

0.053 8.0(11) 1.76(39) 

0.080 8.5(11) 1.96(41) 

0.106 8.9(12) 2.15(44) 

0.132 9.3(12) 2.34(45) 
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APPENDIX G: RESULTS AT THE EXPERIMENTAL QUARK MASSES AT /3 = 6.0 AND 5.85 IN THE 

MS-NDR SCHEME 

Because the MS-NDR scheme is commonly used, we provide here results in the MS-NDR scheme at the individual 
lattice spacings for the _B-parameters and other matrix element ratios at the physical point. These results arc used in 
Sec. VII to obtain our final physical results in that scheme. They are obtained by multiplying the RI/MOM results 
obtained in Sec. VI by U^^^{2GeV)[U^\2GeY)]-\ Here, U^^^ip"^) and U^^p"^) are the matrices which describe 
the running of the matrix elements of interest in the MS-NDR and RI/MOM scheme, respectively. In the case of 
the four-quark operators, we implement them at two loops [57, 58] with Nf — and Aqcd from [59], as in Sec. V. 
The running for the pseudoscalar density is implemented at four loops [60-62], with Nf = and the same value of 
■'^QCD- Our results for all quantities of interest in the MS-NDR scheme at 2 GeV and at the two values of the lattice 
spacing are summarized in Table LV. The systematic errors quoted in these tables are obtained in that same way as 
in Sec. VI. 



TABLE LV: (me + rri) at 2 GeV in the MS scheme as obtained from our simulations at /? = 6.0 and /3 = 5.85. The first error is 
statistical and the second comes from the systematic uncertainty in the determination of Zs- 

13 (m, + mf^{2 GcY) 
6.0 102.(7)(2)MeV 
5.85 107.(5)(8)MeV 



TABLE LVI: Results at the Icaon mass in the MS-NDR scheme at 2 GeV for /3 = 6.0 and 5.85. The first error is statistical and 
the second originates from the systematic uncertainty in the renormalization constants. 





/3 = 6.0 






P = 5.85 




i Bi 


Gi 


Ri 


Bi 




Ri 


1 0.571 (48) (7) 


29.2(26)(4) 


1 


0.541(37) (4) 


27.7(19) (0) 


1 


2 0.679(56)(7) 


-512.(106)(12) 


-17.5(32)(6) 


0.705(71)(33) 


-474. (76) (45) 


-17.1(28)(15) 


3 1.055(77)(44) 


159.(31)(5) 


5.44(97)(9) 


1.143(75)(292) 


154.(18)(13) 


5.56(66)(51) 


4 0.810(41)(11) 


730(105)(8) 


24. 7(32) (3) 


0.781(63)(91) 


624.(84)(29) 


22.5(32)(8) 


5 0.562(39) (12) 


165.(15)(1) 


5.76(73)(3) 


0.517(98) (9) 


140.(28)(18) 


5.0(11)(6) 



TABLE LVII: Results, at /? = 6.0 and 5.85, for the ratios D^^^ of Eq. (16) interpolated to the kaon point and our best estimate 
for ((7r7r)/=2|Q7,8|JS''') in the chiral fimit written as a — » tt"*" matrix element according to Eq. (15), so as to avoid ambiguities 
due to normalization conventions for the two pion state. All results are given in the MS-NDR scheme at 2 GeV. Please see text 
for a discussion of how the central value and systematic error on the matrix elements is obtained. 





D^^^\mIc) [GeV^] 


Dl^^iM],) [GeV^] 




[GeV^] 


(^+\qI'''\k+) 


[GeV^] 


F 


F 


6.0 
5.85 


2.11(20)(1) 
1.79(36) (23) 


9.3(13)(0) 
8.0(11)(4) 


0.211(20)(1)(27) 
0.179(36) (23) (23) 


0.93(13)(0)(12) 
0.80(11)(4)(10) 
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